Lecture 1 



1.1 Introduction 

This course has two aims. 

Aim 1. To explain how Einstein's theory models gravitation as the curvature of space- 
time. This involves setting up some new mathematical machinery, notably tensor calculus. 
The mathematical background will be developed in parallel with the theory. In textbooks, 
it is often done the other way: mathematics first, theory after. But the textbooks do not 
have to be read in linear order. The approach here should make it clearer at each stage 
where we are going. 

Aim 2. To treat the elementary theory of black holes. 

The closest to a course text is Hughston and Tod: An introduction to general relativity 
(LMS student texts 5, CUP). For more advanced reading, R M Wald, General relativity, 
(Chicago) is strongly recommended. [Note: text in red is clickable] 

General relativity is a theory of gravity. It grew from attempts to reconcile 

a) Newton's theory of gravitation 

b) special relativity 

in much the same way (only more quickly) that special relativity grew out of Maxwell's 
theory and Galileo's principle of relativity. As in that case, the reconciliation involved 
rethinking our ideas about space and time — continuing the process begun with SR. 

Special relativity teaches us that we must no longer think of space and time as separate. 
General relativity goes further: it allows space-time to be 'curved', and models gravity as 
the curvature. This leads to some interesting possibilities. For example, the universe might 
be spatially closed like the three-dimensional 'surface' of a sphere in four dimensions. For 
many years, there was some unease about the contrast between such far-reaching conse- 
quences of GR and the rather meagre set of observations through which it was verified — for 
example the fact that the orbit of Mercury rotates by 43" per century faster than it should 
in Newtonian theory and the fact that the positions of stars in the sky are slightly distorted 
by the 'bending of light' in the sun's gravitational field. However this position has changed 
completely in recent years with many new astronomical observations. Perhaps the most 
dramatic of these concern the binary pulsar PSR 1913+16: this consists of a neutron star, 
about 15 miles across, but with a mass of about 1.5 times that of the sun, orbiting another 



1 



star once every 8 hours or so. Here one reverses the Mercury observation: one uses the 
rotation of the orbit to measure the masses, and then calculates the theoretical rate at 
which the orbital period should decrease as the two stars lose energy through gravitational 
radiation. The result, over 15 years, agrees with observation to within 0.5%. 

It is worth keeping in mind this almost unbelievable correspondence between theory 
and observation as we go through what at times will seem some very technical and abstract 
mathematics. In spite of its mathematical abstraction, general relativity is now one of the 
most accurately confirmed theories of all time. An account of the current experimental 
and observational status of general relativity can be seen online in Clifford Wills' article 
in Living Reviews. 

1.2 Newtonian gravity 

The essential content of Newton's theory is contained in two equations. The first is 



where is the gravitational potential p is the matter density, and G = 6.67 x 10 
in SI units (it has dimensions L 3 M _1 T~ 2 ). With appropriate boundary conditions, (1) 
determines the gravitational field generated by a source. The second is 



This determines the force mg on a particle of mass m in a gravitational field with potential 
0. The accuracy of the theory is remarkable: in the solar system, the only detectable dis- 
crepancy between the theoretical and actual motions of the planets is the orbit of Mercury, 
where it amounts to one part in 10 7 . 

The two equations contain the inverse-square law. By integrating over a region bounded 
by a surface S, and containing a total mass M, we obtain Gauss' law from the divergence 
theorem: 



Js 

If the field is spherically symmetric, for example if it is that outside a spherical star, then 
g = F(r)r. By taking S to be a sphere of radius r, we obtain 



V 2 = 47rGp, 



(1) 
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V0. 



(2) 




(3) 



F(r) 



GM 
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which is the inverse square law. The corresponding potential (J) is not unique since we are 
free to add a constant. If we fix this by taking cf) = at infinity, then 



E 



\mv 2 + mcj) 



= \mv 2 — 



GMm 



r 



is the total energy of a particle at of mass m and speed v. 

One can see even in these very simple formulas some important issues. 

Black holes. For the particle to escape to infinity, v must exceed the escape velocity 



(since E is conserved and v 2 must remain non-negative as r — > oo). This is maximal at 
the star's surface, where r takes its lower possible value in the region outside the star. 
What if at this point we have v = c, the velocity of light? This will be the case if the 
radius of the star is R = 2GM/c 2 , the so-called Schwarzs child radius. Then nothing can 
escape from the surface. But what if there is a mirror on the surface and we shine light 
down from infinity? It will be reflected at the surface and follow the same path back out 
again. Since orbits are reversible in Newtonian gravity, it will be reflected back to infinity. 
Clearly Newton does not provide a consistent picture of such a 'black hole'. One might 
not think this matters: for a star of the mass of the sun, the Schwarzschild radius is 3 
km, which would make the star an apparently impossibly dense object with a density 10 16 
times that of water. However, such stars exist not much more massive than the sun. If 
one increases the mass to that of a galaxy, then the critical density goes down to that of 
air. So the issue must be addressed as a potentially real physical problem, and not simply 
as an absurd extrapolation. 

Linearity. The equations of Newtonian gravity are linear, if you superimpose two mass 
distributions, then you need only add the corresponding gravitational fields. However, as 
bodies interact gravitationally, energy is transferred from their gravitational fields to the 
bodies themselves, and vice versa. Thus gravitational fields carry energy, and therefore by 
Einstein's principle E = mc 2 , they have mass: that is gravitational fields should themselves 
generate gravitational fields. Thus a relativistic theory of gravity should be described by 
nonlinear equations. 



v = sJlGMj* 



r 
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1.3 Gravity and SR 



(1) and (2) are very like the basic equations of electrostatics. This suggests that we 
should try to construct a relativistic theory of gravity in the same way as you did for 
electromagnetism. But nobody has succeeded in doing this in a way that works — i.e. 
that does not give incorrect predictions; and there is any case problems of principle about 
applying special relativity in a gravitational field. 

Gravity and light.. The first problem is that light is affected by gravity, and this is 
incompatible with the basic notion in special relativity that light should travel in straight 
lines at the same constant speed in all inertial frames. 

How do we know that light is affected by gravity? First because it is observed to be. 
For example, the path followed by light reaching the earth from a star in the direction of 
the sun is bent by the sun's gravitation field, as was first observed by Eddington's in 1919: 
he showed that, relative to their normal positions in the night sky, the positions of stars 
appeared to be displaced away from the centre of the sun when observed during a total 
eclipse. (The history of Eddington's observation is not quite as straightforward as it is 
sometimes presented. See, for example, Peter Coles' article Einstein, Eddington and the 
1919 eclipse. More recent and dramatic evidence for the bending of light can be seen in 
the photographs taken by the Hubble telescope.) 



Second, there is Bondi's thought experiment: if photons were not affected by gravity, 
then one could in principle build a perpetual motion machine. Bondi imagined a machine 
consisting of a series of buckets attached to a conveyor belt. Each contains a single atom, 




a Apparent 



Sun 



Figure 1 : Bending of light 
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with those on the right in an excited state and those on the left in a lower energy state. 
As these reach the bottom of the belt, the excited atoms emit light which is focused by 
two curved mirrors onto the atom at the top of the belt; the one at the bottom falls into 
the lower state and the one at the top is excited. Since E = mc 2 , those on the right, 
which have more energy, should be heavier. The force of gravity should therefore keep 
the belt rotating. The way out is that photons lose energy as they climb up through the 




Figure 2: Bondi's perpetuum mobile 

gravitational field. Since E = hu, this means that they are red-shifted. Such a shift was 
measured directly by Pound and Rebka in 1959 by a remarkable experiment: over the 75ft 
height of the tower of the Jefferson building at Harvard, the red-shift is about 3 parts in 
10 14 . 

The Pound-Rebka result is incompatible with SR, as illustrated in the space-time diagram. 
The blue lines are the top and bottom of the tower and the red lines are the worldlines 
of photons travelling up the tower. Since the top and bottom of the tower are at rest 
relative to each other, their worldlines in SR are parallel, which forces At = At'. So in a 
special-relativistic theory of gravity, there cannot be any gravitational red-shift. 

Where are the inertial frames?. In special relativity, there is a preferred class of 
'non-accelerating frames', called inertial frames. These are precisely the frames in which 
Newton's first law holds: 

In the absence of forces, particles move in straight lines at constant speed. 
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Figure 3: Pound and Rebka's measurement is incompatible with special relativity 

How can we tell if we are in an inertial frame? At first glance the answer looks easy: we 
simply look to see if free particles travel in straight lines at constant speed relative to the 
frame. The problem is that gravity affects all matter equally, so there are no completely 
free particles. In the equation of motion, 



the two m's cancel, so all particles accelerate in the same way. This does not happen with 
other forces. In electromagnetism, there are (charged) particles that are affected by an 
electric field and others (uncharged) that are not; but in gravitation theory, there are no 
'neutral' particles which we can think of as free of all forces. 

The m on the LHS is inertial mass; that on the RHS is gravitational mass (the analogue 
of charge in electromagnetism). The exact equality of the two is the (weak) principle of 
equivalence. Galileo tested it by comparing the periods of pendula with weights made out 
of different substances. A celebrated 19 th experiment by Eotvos checked very accurately 
(to one part in 10 9 ) the equality of the two ms in the gravitational force on a mass and the 
centrifugal force due to the Earth's rotation: these are the two components in the apparent 
gravitational field. More recent experiments — including lunar ranging measurements — 
have improvement the accuracy to one part in 10 13 ; and a planned space experiment 
(STEP) will test it to one part in 10 18 . 

We have good reason, therefore, to believe that the equality is exact. Einstein extrap- 
olated it to the (strong) equivalence principle: there is no observable distinction between 
the effects of gravity and acceleration. In a closed box, there is no physical experiment 



mr = — mV0 



(4) 
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Figure 4: Relative acceleration in free fall 



within the box that will reveal whether(i) the box is at rest on the earth's surface, or (ii) 
the box is accelerating at lg in otherwise empty space. In both situations, those inside the 
box 'feel' a normal terrestrial gravitational field. 

In a gravitational field, it is impossible to identify the global inertial frames of special 
relativity. But we can pick out local inertial frames in which gravity is turned off: these 
are the frames in free-fall. But if we can only work in frames in which gravity is turned off, 
then how can we observe gravity? The answer is: by shifting your point of view. Gravity 
manifests not as a 'force', but rather as the relative acceleration of nearby local inertial 
frames. The key difference between what happens near the earth and what happens in 
empty space is that (i) near the earth, two nearby 'boxes' in free fall (e.g. space stations) 
have a small relative acceleration because they are falling on paths that converge at the 
centre of the earth, while (ii) in empty space, they have no relative acceleration. 

Starting point. Special relativity holds over short distances and times in 
frames in free fall. Gravity is not a local force field, but shows up in the small 
relative acceleration between local inertial frames. In the presence of gravity, 
the transformation between local 'inertial' coordinates is not exactly linear. 

The idea of curvature comes in here: the analogy is with map making. If one makes maps 
of the earth's surface by projecting onto a tangent plane from the centre of the earth, then 
overlapping maps will be slightly distorted relative to each other because of the curvature 
of the earth. 
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Lecture 2 



2.1 Inertial coordinates 



SR describes the relationship between physical observations made by different non-accelerating 
observers. Each observer labels events in space-time by four coordinates t, x, y, z. The co- 
ordinate systems of two observers are related by a Lorentz transformation 



(t\ 
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where T is a column vector (which shifts the origin of the coordinates) and 



(5) 
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(6) 
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with L° > 



detL = 1, and L^ 1 = gL l g, where g is the diagonal matrix with diagonal 
entries 1,-1,-1,-1, (i.e. L is a proper orthochronous Lorentz transformation matrix). 
For example for a boost along the x-axis, 
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where 7 = 1/v^l — u 2 . Throughout this course, we shall take c = 1. 
We call t, x, y, z an inertial coordinate system or an inertial frame. 
To make the transition to GR, we need some new terminology and notation. 

Inertial coordinates. These will be written a; = t, x 1 = x, x 2 = y, x 3 = 
important to keep track of the position of the indices: a lot of information will 
by making a distinction between upper and lower indices). We can then write 
compact form 



z. (It is 
be stored 
(5) in the 



x 



L\x b + T a 



6=0 



0,1,2,3). 



(When you write the entries L a b in L, make it clear that the upper index comes first: a 
labels the rows in L and b the columns). We note that 



_ dx a 
L b ~ dx» 

and that 

dec ^ 

(L ~ ln = dx~» • 

The summation and range convention. When an index is repeated in an expression 
(a dummy index), a sum over 0,1,2,3 is implied. An index that is not summed is a free 
index: any equation is understood to hold for all possible values of its free indices. To 
apply the conventions consistently, an index must never appear more than twice in any 
term in an expression, once as an upper index and once as a lower index. 

The metric and the Kronecker delta. We define the quantities g a b, g ab and 5% by 



1 a = b = 
</„/, =g ab = { -1 a = b^0 5 a b 
otherwise 



1 a = b 
otherwise 



(Note: later on, the 'metric coefficients' g a b and g ab will be more general functions of the 
space-time coordinates, but the Kronecker delta will always be defined in this way.) 

Examples. (8) becomes 

x a = L\x b + T a . (9) 

(repetition of b implies summation over 0,1,2,3, while the range convention means that 
the equation is understood to hold as the free index a runs over the values 0, 1, 2, 3). If two 
events have coordinates x a and y a in the first system and x a and y a in the second system, 
then 

x a - y a = L\(x b - y b ) = L\x b - L\y b . (10) 

(This illustrates that one must take care about what is meant by a 'term in an expression'. 
In principle, you should multiply out all the brackets before applying the rule — otherwise 
the three- fold repetition of b in the middle expression could cause confusion). The Lorentz 
condition L l gL = g becomes 

d dx c dx d 

L a L b9cd = 9"lQ^Q^b = 9ab ■ 
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(Note that it does not matter in which order one writes the Ls and gs so long as the indices 
are 'wired up' correctly. In this equation a, b are free, while c, d are dummy indices — like 
dummy variables in an integral. The sum over c is the sum in the matrix product tfg, 
while the sum over d is the sum in the matrix product gV). Similarly, L~ x g~ 1 {V'')~ 1 = g~ x 
becomes 

If one combines two coordinate transformations 

x a = K a b x\ x a = L\x b + T a (12) 

then the result is 

x a = K\L b c x c + K\T b . (13) 

(To avoid ambiguity, one changes the dummy index in the second equation before making 
the substitution. It is then clear that there are two sums, over b = 0,1,2,3 and over 
c = 0, 1, 2, 3. If you did not do this, then you would end up with the ambiguous expression 
K a b L b b , which could mean Y^l=o ^b^b)- A a final illustration, we note the identity 

9ab9 bc = S c a - (14) 
(This is the matrix identity g 2 = 1. Here c, a are the free indices and b is a dummy index). 



2.2 4- vectors 



A 4- vector has four components V°, V 1 , V 2 , V 3 , which transform under change of coordi- 
nates by 

/V°\ /V°\ 



V 1 
V 2 



L 



V 1 
V 2 

\v 3 J 



(15) 



That is V a = L a b V b . Note that L a b = dx a /dx b . The following definition will generalize 
easily to GR. 

Definition 1 

A 4- vector is an object with components V a which transform by 

dx a ~ 



ya 



dx b 



V° 



under change of inertial coordinates. 
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Example. The 4- velocity: if x a = x a (s) is the worldline of a particle, parametrized by 
proper time s, then the 4- velocity has components V a = dx a /ds. Under coordinate change 



ya 



dx a dx a dx b 



ds dx b ds ' 

so the 4-vector transformation rule is a consequence of the chain rule. 



(16) 



2.3 Tensors in Minkowski space 

Two other objects have similar transformation rules. 
(1) The electromagnetic (EM) field 
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\F 30 


F 31 


F 32 


F 33 J 



transforms by F = LFL f . That is 



nab 



L\L\F 



cd 



dx a dx b p cd 
dx c dx d 



(17) 



(2) The gradient covector of a function f(x a ) has components d a f, where d a = d/dx a . 
These transform by the chain rule 



Oaf 



dx ~ 



dx a 



(19) 



Note that it is dx/dx on the RHS, not dx/dx, so this is not the 4-vector transformation 
rule. Hence the term 'covector'. 

Definition 2 

A tensor of type (p, q) is an object that assigns a set of components T a --- b c . ^ (p upper 
indices, q lower indices) to each inertial coordinate system, with the transformation rule 
under change of inertial coordinates 



Ta...b 
c.d 



dx a dx b dx h 



dx e dxf dx c dx d 



J- h...k- 
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A tensor can be defined at a single event, or along a curve, or on the whole of space-time, 
in which case the components are functions of the coordinates and we call T a tensor field. 
If q — then there are only upper indices and the tensor is said to be contravariant; if 
q = 0, then there are only lower indices and the tensor is said to be covariant. 

Examples. (1) A 4-vector V a is a tensor of type (1, 0). Also called vector or contravariant 
vector. 

(2) The gradient covector d a f is a tensor of type (0,1). A tensor a a of type (0,1) is 
generally called a covector or covariant vector. 

(3) The Kronecker delta is a tensor of type (1, 1) since 



(4) The contravariant metric has components g ab and is a tensor of type (2,0) (by 11). 
The covariant metric has components g a t, and is a tensor of type (0, 2) 

Both the Kronecker delta and the metric are special in that they have the same com- 
ponents in every frame. For a general tensor, the components in different frames are not 
the same. 



dx a dx d 
d dS^dx E 



dx a dx' 



dx c dx b 




(20) 
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Lecture 3 



3.1 Operations on tensors 

Addition. S, T of the same type: S + T has components 

oa...b i rpa...b 

J c...d + 1 c...d- 

Scalars. A scalar at an event is simply a number. A scalar field is a function on space- 
time. The value of a scalar is unchanged by coordinate transformations. We can multiply 
a tensor T by a scalar / to get a tensor of the same type with components /T a " b c ...d. 

The two operations make the space of tensors of type (p, q) into a vector space of dimension 
A p+q . 

Tensor product. It S, T are tensors of types (p,q), {r,s), respectively, then the tensor 
product is the tensor of type (p + r, q + s) with components S a '" h ' c ...dT e '"^ g ...h- It is denoted 
by ST or S ® T. 

Differentiation. If T is a tensor field of type (p, q), then VT is the tensor of type (p, q+1) 
with components 

be be ^ 

V a T "' C d...e = 9 a T "' C d...e, d a 



Under change of inertial coordinates, 



d a T b - c d e = ^-d, 

Q x a > 



dx a ' 



dx h dx r ^ q 
dx q dx d 



8X P dx b q rpq... 

dx a dx q dx d p 

which is the correct transformation rule. Note that we are still working in the context of 
special relativity: the calculation only works because dx/dx is constant. We shall have to 
work harder to define differentiation in curved space-time. 

Contraction. If T is of type (p + 1, q + 1), then we can form a tensor S of type (p, q) by 
contracting on the first upper index and first lower index of T: 

ob...c _ rpab...c 

O ,.../ - J- ae...f 
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(sum over a). Under change of coordinates 



e...f 



Tab...c 
ae...f 

dx a dx b 



dx c dx s dx 



dx p dx q dx T dx a dx e dx-f 



st...u 



dx b 
dx q 



since 



dx° dx 
dx r dx e 

dx a dx s 



dx? 



5 s . 

u p . 



t...u 



dx p dx a 

One can also contract on other pairs of indices (one upper and one lower). 

Raising and lowering. If a is a covector and U a = g ab ab, then U is a 4-vector (tensor 
multiplication plus contraction). We write a a for U a and call the operation 'raising the 
index': raising the index changes the signs of the 1,2,3 components, but leaves the first 
component unchanged. The reverse operation is 'lowering the index': V a = g a bV b . One 
similarly lowers and raises indices on tensors by taking the tensor product with the covari- 
ant or contravariant metric and contracting, e.g. T a b = gi, c T ac . Warning: keep track of the 
order of the upper and lower indices: T\ and T b a are generally distinct, so don't write T b a . 

Examples. (1) If / is scalar field, then V a /, where 

(VV) = (dj, -d X f, -dyf, -dj) 

is a 4-vector field: it the 'gradient 4-vector'. 

(2) If U and V are 4-vectors, then 

g(u, V) = g ab U a V b = u a v a = u a v a . 

(3) Raising one index on g a & or lowering one index on g ab gives the Kronecker delta since 



g gbc - o c . 



(4) Suppose that (S a ) 
respective components 



(S a T b 



[1, 0, 0, 0) and (T a ) = (1, 1, 0, 0). Then S ® T and T <g> S have 



(\ 






1 

















0/ 



(T a S b 



(I 

1 






















0/ 



Note that S ( 
S ® T and T 



T ^ T ® S, but when written as matrices, as above, the components of 
) S are related by transposition. 
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3.2 The energy- momentum tensor 



Consider a cloud of particles ('dust'), in which the velocities of the individual particles 
vary smoothly from event to event and time to time. There is one worldline through each 
event and the 4-velocity is a 4-vector field. There are no external forces or interactions, so 
each particle moves in a straight line at constant speed. 



Problem What is the energy density seen by an observer with 4-velocity VI The answer 
depends on V because 

(i) the energy of each individual particle depends on its velocity relative to the 
observer; and 

(ii) moving volumes appear to contract. 
Definition 3 

The rest density p is a scalar. It is defined at an event A to be the rest-mass per unit 
volume measured in a frame in which the particles at A are at rest. If there are n particles 
per unit volume in this frame and each has rest mass m, then p = nm. 




Space 



Space-time 



Figure 5: A 'dust' cloud 
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Consider the particles which occupy a unit volume at an event A in the rest frame of the 
particles at A. Suppose that in this frame the observer is moving along the negative x-axis 
with velocity v. To the observer, each particle at A appears to have velocity (v, 0, 0) and 
to have energy 1717(1)) = mj \J\ — v 2 . The particles appear to occupy a volume 1/7(17) = 
\J\ — v 2 . Therefore the observer measures the energy density to be 7(u) 2 p. 

Definition 4 

The energy-momentum tensor of the dust cloud is the tensor field with components T ab = 
pU a U b . It is a tensor of type (2, 0) (scalar times the tensor product of two 4- vectors). 




Figure 6: The transformation of density 

Proposition 1 

The energy density measured by an observer with 4-velocity V is p Q b s = T ab V a V b . 
Proof. In the rest frame of the observer, 

(V a ) = (1,0,0,0) (Z7 a )= 7 (v)(l,v,0,0). 

Therefore T ab V a V b = p(U a V a ) 2 = p-f(v) 2 . 
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3.3 The continuity equation 

In the rest frame of the particle P, the nearby particles all have small velocities and so the 
non-relativistic continuity equation holds 

^+pdivu = 0. 

This expresses the conservation of mass and energy. 

Proposition 2 
W a T ab = 

Proof. We want to show that this holds at an event A. It is a tensor equation, and so it 
enough to show that it holds in one frame. Now 

v a T ab = v a ( P u a u b ) = (u a v aP + P (y a u a ))u b + pu a v a ii b . 

Each particle moves in a straight line at constant speed. Therefore 

Q= dU^ = dx^dU^ = ubVbU ^ 
ds ds dx b 



In the rest frame of the particle at A 

U = 



(1,0) at A 
7(u)(l,u) near A, u — |u| << 1 



Therefore, at the event A, 



u a v aP = d tP 

V a U a = d tl {u) + d x (j(u) Ul ) + d y (j{u)u 2 ) + d z (j(u)u 3 ) , 

where d t = d/dt and so on. But since ^(u) = 1/ ' \J\ — u 2 , we have 

dt"f(u) = ^(ufudtu = 

at A, where u = 0. The other partial derivatives of 7(w) vanish similarly. Therefore 
V a U a = div u at A and the result follows from the continuity equation. ■ 
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Lecture 4 



4.1 Curved space-time 

In Minkowski space, the components of the metric tensor are constant. In general relativity, 
they are functions of the space-time coordinates. The gravitational field is encoded in the 
second derivatives of the metric, which measure the curvature of space-time. Our starting 
point is the following set of guiding principles. 

GR1: SR holds over small distances and short times in frames in free fall, that 
is, in local inertial frames. In such frames we can set up local inertial coordinates 
and as in Minkowski space. 

GR2: Gravity appears as the relative acceleration of nearby local inertial 
frames. 

We can only use inertial coordinates in a small neighbourhood of each event. If we want 
to describe what is happening in a larger region of space-time, we must use a general 
coordinate system x a (a = 0, 1, 2, 3) to label events. 

GR3: Principle of general covariance: the theory should be formulated in a way 
that does not give any special role to a particular coordinate system. 

This generalizes the principle of relativity, which says the same for inertial coordinate 
systems in SR. 

4.2 The metric 

Let A be an event and let O be an observer in free fall at A. To O, SR must hold near 
A. Suppose that A sets up local inertial coordinates t, x, y, z in small region of space-time 
around A, with A as origin. Since O thinks SR holds locally, we have the following. 

(1) The worldlines of free particles that pass through A satisfy 



d 2 t d 2 x d 2 y 



d 2 ^ 
d^ 







ds 2 ds 2 ds 2 



sAt — x — y — z — 0, where s is proper time. 
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(2) Let B be an event near A with coordinates dt, dx, dy, dz, and put 

ds 2 = dt 2 _ dx 2 _ dy 2 _ dz 2_ 

Then, provided that we ignore third order terms, ds 2 has the same interpretation as in SR. 
Therefore the following hold. 

Timelike: If ds 2 > 0, then ds is the time from A to B on a clock in free fall 
from A to B. 

Null: If ds 2 = 0, then A and B lie on the worldline of a photon. 

Spacelike: If ds 2 < 0, then ds 2 = — D 2 , where D is the distance from A to B 
measured in a frame in free fall in which A and B are simultaneous. 

The only difference from SR is that in SR we would say that (1) held everywhere and that 
(2) was also true for widely separated events. 




Figure 7: The displacement from A to B in the three cases 

The coordinates t, x, y, z can only be used at A. Let us introduce a general coordinate 
system x a on space-time, and translate these statements into the general system. Near A, 
t, x, y, z are some general functions of the x a s, so 

dt 

dt = — — dx a + second order terms in dx 
ox a 

and so on. 
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Therefore at A, ignoring third order terms in the dx a s, 

ds 2 = g ab dx a dx\ 

where, in a generalization of our previous notation, the coefficients g a b = gba are functions 
of the coordinates constructed from the partial derivatives of the local inertial coordinates 
at an event with respect to the x a s. They are called the metric coefficients. For example, 
at A, 

dt dt dx dx dy dy dz dz 
901 dx° dx 1 dx° dx 1 dx° dx 1 dx° dx 1 
We can make a similar transformation to local inertial coordinate systems near other events. 
Then ds 2 = g a bdx a dx b everywhere, but the g a bS now vary from event to event. 
If we change from x a to a new coordinate system x a , then 



^ = ^d/=^gg) dx°dx b . 



In the new coordinate system, the metric coefficients are 

dx c dx d 
9ah = 9cd W a dx~ b 

or in matrix notation 

g = M l gM where M = 



dx a 
dx b 

At any one event, we can reduce g to the the diagonal matrix with diagonal entries 
1,-1,— 1,-1 by transforming to local interial coordinates at that event (but in general we 
cannot do this at all events simulataneously, as we shall see). Therefore the matrix g has 
one positive and three negative eigenvalues — which is usually expressed by saying that g 

has signature H . 

To summarize, in an arbitrary coordinates system, if dx a is the coordinate separation 
between two nearby events A and B, then, to the second order in da; 

ds 2 = g ab dx a dx b 

where the metric coefficients are evaluated at A and ds has the interpretation in (dsl)- 
(ds3). 

The metric cooemcients g ab have the following properties 
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MCI 



they are smooth functions of the coordinates x a ; 
they are symmetric g ab = g ba 

they have signature H at every event 

they transform under general coordinate transformations by 



MC2 



MC3 



MC4 



dx c dx d 



9ab — 9cd 



dx a dx b 



Example. Suppose that a; = t, x l = r, x 2 = 8, x 3 = (f) and 



ds 2 = dt 2 - dr 2 - r 2 d9 2 - r 2 sin 2 Odcf 



2 



Then we can reduce ds 2 to the form dt 2 — dx 2 — dy 2 — dz 2 by the coordinate change 
t = t, x = rsin^cos^, y = rsin#sin0, z = rcos8. So this is just the metric of SR in a 
funny coordinate system (spherical polars). For a general metric, we cannot reduce to the 
Minkowski space form by a coordinate transformation. However, we can do it up to the 
second order at one event: the following proposition will allow us to recover local inertial 
coordinates from a general set of metric coefficients satisfying (MC1)-(MC4). 

Proposition 3 

Let g a f, be a general set of metric coefficients such that (MC1)-(MC4) hold and let A be 
the event x a = 0. Then there exists a coordinate system x a such that x a = and d c g a b = 
at A. 

Proof. Define new coordinates x a by c , where the r a bc s are constants 



such that r a bc = T a cb . Put G a b = g a bi^) and H cab = d c g a b(0). Then, by Taylor's theorem, 

g ab = G ab + x c H cab + 0{2), 
where '0(2)' denotes quadratic and higher order terms in the x a s. Hence 

dx° dx d 



g a b — 9cd 



dx a dx b 



(G cd + x k H kcd )(S c a - T c ae x e )(S d b - + 0(2) 

Gab + x c {H cab — r 

abc r?>ac ) + 0(2), 
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where T abc = G ad T d bc ; we have used x a = x a + 0(2), as well as changing the names of the 
dummy indices. We want to choose Y abc = T acb so that 

H C ab Tabc Tfoac • 

By permuting the indices, we then also have 

H bca r ^ cab ~\~ I" \ cb H abc Tfeca ^ cba ■ 

By adding the first two of these and subtracting the third, we obtain 

Tabc 2^^ ca ^ ^bca H abc ^ i 

and hence that 

ra _ XriadiTT \ zj IT \ 

be — 2^ \ n cdb + tl hcd — ti dbc ) , 

where G ab Gb c = 5% (that is, (G ab ) is the inverse of the matrix G ab . 
Conversely, if we define Y a bc in this way, then we have 

H cab — ^abc — F bac = H cab — \(H cab + H bca — H abc + H cba + H acb — H bac ) = 0, 

since H abc = H acb . ■ 
Note that 

TV = \g ad {d c9db + d b9dc - d d9bc ) . 

evaluated at x a = 0, where the g ab s are the inverse (or contravariant) metric coefficients, 
defined by g ab g bc = 5%. The quantities T bc are called the Christoffel symbols. We shall meet 
them again in the definition of the Levi-Civita connection. 
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Lecture 5 



5.1 Existence of local inert ial coordinates 

We are trying to extract local inertial coordinates from a general set of metric coefficients 
g a b satisfying (MC1)-(MC4). We have shown that by coordinate change, we can make 
d a gbc = at the event x a = 0, that is 

9ab(x) = g ab {0) + O{2). 

Proposition 4 

Let A be an event. Suppose that we have two coordinate systems x a and x a such that 
at A, x a = x a = and d a gb c = d a 9bc = 0. Then there exist constants M\ such that 
x a = M\x b + 0(3) ('0(3)' means third order terms in x). 

This says that the transformation is linear at A up to the second order in x: that is, the 
Taylor expansion about A of x a in powers x a has no second order terms. 

Proof. We have to show that at A, d 2 x a / dx b dx c = 0. Now at all events, 

dx c dx d 



9ab 9cd ~„ rs-f. ■ 

OX a OX 

Therefore 

_ dx? dx c dx d d 2 x c dx d dx c d 2 x d 

e 9ab Q~ e f^ cd g^aQ^b cd Qx a 8x e 8x h dx a dx b 8x e 

(note that the second two terms on the right-hand side differ by the interchange of a and 
b). Put 

dx c d 2 x d 



J abe 9cd r>~ 



dx a dx b dx e 



Then L abe = L aeb and (21) gives 



L ba e + L abe = 
L e ba + L bea = 
L a eb + L eab = . 

By adding the first and third, and subtracting the second, we obtain L abe = 0. Hence 

d 2 x q dx a „ n dx c d 2 x d dx a „ nr 

9 pq 9cd^-^j^- = ^9 Pq L abe = , 



dx b dx e dx p dx a dx b dx e dx p 
which completes the proof. 
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Proposition 5 

Existence of local inertial coordinates for a general metric. Let g a b{x) be an set of metric 
coefficients satisfying (MC1)-(MC4) and let A be an event. Then there exists a coordinate 
system x a such that x a = at A and 



(9ab(x)) 



(I 








o \ 





-f 














-f 





\o 








-v 



+ 0(2). 



This system is unique up to coordinate transformations x a 
for some Lorentz transformation matrix L = (L a b ). 



x a where x a 



T a ~b 
L, b X 



+ 0(3) 



(I 








o \ 





-1 














-1 













-1/ 



Proof. Choose an initial coordinate system such that d c g ab = and x a = at A. Since 
g a b has signature H we can find a matrix M = (M a b ) such that 



M l gM = (M a c g ab (0)M b d ) 



Now make the linear coordinate change x a h- > M a b x b . The uniqueness statement follows 
from the previous proposition. ■ 

The metric coefficients g a b will describe a general gravitational field. The coordinates at 
A in the last proposition are interpreted as local inertial coordinates of an observer in free 
fall at A. For special metrics we can reduce g ab to the diagonal form diag (1,-1,-1,-1) 
everywhere. We shall see that that is the case in which gravitational field vanishes, but for 
a general metric, such a coordinate transformation does not exist. To summarize: 

(1) A gravitational field is described by a general set of metric coefficients 
satisfying (MC1)-(MC4), which encode the temporal and spatial separation of 
nearby events according to (dsl)-(ds3). 

(2) The local inertial coordinates set up by an observer in free fall at an event 
A are the coordinates x a such that at A, x a = and 



(9, 



ab) 



(I 








o \ 





-1 














-1 













-v 



0(2). 
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In local inertial coordinates, SR holds over small times and distances. It is conventional 
to write the metric in general coordinates in the form 

ds 2 = g ab dx a dx b . 

For example, for Minkowski space in spherical polar coordinates 

ds 2 = df 2 _ dr 2 _ r 2 dQ 2 _ r 2 gin 2 g & , 2 _ 



5.2 Particle motion 

In a local inertial coordinate system at an event A, d c g a b = at A and the worldlines of 
free particles (particles in free fall) satisfy 

x a = where ■ = — , s — proper time 
ds 

at the event A. What is the corresponding differential equation in a general coordinate 
system? 

Our strategy is find a Lagrangian and use a result from classical mechanics about the 
transformation of Lagrange's equations under change of coordinates. 

Invariance of Lagrange's equations. In classical mechanics, the equations of motion 
of a system with time-independent Lagrangian L(q a ,q a ) are Lagrange's equations, 

d { dL\ dL 

0, 



dt \dq a J dq a 

where the q a s are generalized coordinates. The equations in a new coordinate system q a 
can be found by substituting 

Qa = q a {q), q a = ~^rqb 

into L and writing down Lagrange's equations in the new coordinates. This is the sense in 
which Lagrange's equations are invariant under coordinate transformations. 
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We shall take this result out of its original physical context and apply it as a statement 
about the way in which Lagrange's equations change under change of dependent variables. 
We shall put the space-time coordinates x a in the role of the q a s and the proper time s in 
the role of time in classical mechanics. We are led to the following. 

The geodesic hypothesis.. The worldlines of particles in free fall parametrized by proper 
time s satisfy the geodesic equations 

d (SL\_BL =0 
ds \dx a J dx a 

where • = d/ds and L = \g ab x a x b . Solution curves are called geodesies. 

Justification (1). In local inertial coordinates at an event, these reduce to 
x a = at the event. 

Justification (2). They are invariant: they take the same form in every coordi- 
nate system. 

To understand (1), we write out the geodesic equations explicitly. We obtain 

^(g ab x b )-l(d a g bc )x b x c = 0. 

Hence 

g a bX b + \x h x c (2d c g ab - d a g bc ) = 
=> x a + \x b x c g ad {d b g dc + d c g bd - d d g bc ) = 

— I* x ~~\~ r^ c x x —— o 

where 

T tc = \g ad {db9dc + d c g bd - d d g bc ). 

These are the Christoffel symbols or connection coefficients. We shall see a lot of them. 
Since the Christoffel symbols vanish at an event A in local inertial coordinates at A, the 
equations reduce to x a = in these coordinates at A. 

The second statement in justification follows from the general principles of Lagrangian 
mechanics and from the invariance of L: from (MC3), 

_ ^ b _ _ dx^dx^_. c . d _ . c . d 

g a b% X 9ab r> „ n ^ 9cd-E % ■ 

OX c OX d 
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Putting these two facts together, we have that the geodesic equations give the correct 
equations of motion at any chosen event in one particular coordinate system; and also that 
if they hold in one system, then they hold in every system. Hence they give the correct 
equations of motion everywhere in every coordinate system 

It follows from the geodesic equations that L is constant. In fact, on the worldline of a 
particle in free fall, ds 2 = g a bdx a dx b , by definition, so L = g a bX a x b = 1. 

Example. For the Minkowski space metric written in spherical polar coordinates, 

L = \(i 2 - r 2 - r 2 6 2 - r 2 sin 2 6<p 2 ) . 

The geodesic equations are 

i = 6 + 2r~ 1 f6 -sm6cos6<p 2 = 

r-rQ 2 -r sin 2 Ocf 2 = 0, <j> + 2r~ 1 f</> + 2cot ^0 = 0. 

We can read off from these that, for example, r^ 3 = 1/r (we are labelling the coordinates 
so that x° = t, x 1 = r, x 2 = 6, x 3 = 0). 

We can extend the discussion to photons, where the same argument leads to the pho- 
ton postulate that the worldline of a photon is also given by the geodesic equations with 
g a bX a x b = 0. The parameter s does not here have the interpretation of time: it is called 
an affine parameter and can be replaced by any linear function of s. Geodesies with 
g a bX a x b > are said to be timelike] those with g a bX a x b = are said to be null. 
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Lecture 6 



6.1 Transformation of the Christoffel symbols 

The Christoffel symbols are defined by 

T lc = \g ad {d b g dc + d c g ba - d a g bc ). 
If we make a coordinate change, then 

r& c = \g ad {db9dc + d c g ba - d a g bc ). 
How are T bc and T bc related? We know that the geodesic equations 

x a + T a bc x b x c = (22) 

transform to 

l a + f a bc x b x c = 

because the Lagrangian they are derived from is invariant. Substitute 

5* = —x d 
dx d 

into the second equation to get 

dx a ,. d d 2 x a , d . e ~ dx e dx f . b . c 

U = ~,X ~\- — -77- X X ~\- 1 p f~^Z 

ox d ox a ox e J ox b ox c 



~ j dx e dx? d 2 x d 



°f dx h dx c dx h dx c 



x b x c , (23) 



with the second line following form the first by multiplying by dx p / dx a and summing over 
a. Hence since T bc = T® b , and since (23) and (22) are equivalent for all choices of free 
particle worldline, 

dx a ~ j dx e dx? dx a d 2 x d 
bc ~ dx d ef dx b dx c dx d dx b dx c 
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6.2 Manifolds 



We now have half of general relativity. 

(1) The gravitational field is encoded in the metric coefficients g a b. 

(2) Free particles move on timelike geodesies; photons move on null geodesies. 

The picture is very like SR except that the light cones vary from event to event. So 
we know how gravity affects matter. We shall leave the other half — how matter generates 
gravity — until later. First we need some mathematical ideas to make precise our model of 
space-time. 

What sort of object is space-time? We want to draw an analogy between the metric and 
the first fundamental form of a surface. A surface in space can be described parametrically 
by r = r(u, v), where u, v are parameters. The distance ds between the nearby points 
(u, v) and (u + du, v + dv) is given by 

ds 2 = r u .r u du 2 + 2r u .v v dudv + v v .v v dv 2 
= Edu 2 + 2Fdudv + Gdv 2 

This is the first fundamental form (never mind about the second). Like ds 2 in space-time, 
it is a quadratic form in (du, dv), which measures the separation between two nearby points 
on the surface. The coefficients E, F, G are functions of the 'coordinates' u, v, like the 
metric coefficients g a j, in space-time. We note two points. 

(f) In general, the metric cannot be reduced to the flat form du 2 + dt> 2 by 
changing the parameters (this is only possible if the surface has no intrinsic or 
Gaussian curvature). 

(2) The surface may have nontrivial topology, in which case the parametrization 
cannot be used over the whole surface. In GR, similarly, we must allow for 
space-time to have a nontrivial topology: this plays an important part in the 
model space-times used in cosmology. 

We want space-time to be a 'four-dimensional surface', but we don't want to have to think 
of it as embedded in a higher-dimensional flat space. Locally, the surface looks like a piece 
of 'curved Minkowski space' just as the surface locally looks a 'curved piece of the plane'. 
The appropriate mathematical object is the following. 
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Definition 5 

A four-dimensional manifold is 

(a) a connected Hausdorff topological space M, together with 

(b) a collection of charts or coordinate patches (U, x a ), where U C M is an open 
set and the x a s are functions x a : U — > R, a — 0, 1, 2, 3, such that the map 

x : U -> K 4 : in h (x° (m) , x 1 (m) , x 2 (m) , x 3 (m)) 

is a homeomorphism from £7 to an open subset Kcl 4 . 

Two conditions must hold (i) every point (='evenf) of M must lie in a coordinate patch; 
and (ii) if (U, x a ) and (U, x a ) are charts, then the x a s can be expressed as functions of the 
x a s on the intersection. We require that 

( ™0 ™1 ™2 ™3\ / 1 2 ™3\ 

should be smooth (i.e. infinitely differentiable) and one-to-one, with 

dx a 
dx b 

The same definition works, of course, in other dimensions. 
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6.3 Vectors and tensors 

The various physical objects in space-time are represented by scalars (functions on space- 
time) or by scalars or tensors, which are objects which transform in simple ways under 
change of coordinates. The definition is the same as in SR, except that the coordinate 
changes are now general. 

Definition 6 

A tensor of type (p, q) is an object that assigns a set of components T a - b c d (p upper 
indices, q lower indices) to each local coordinate system, with the transformation rule 
under change of coordinates 

rpa...b d xCL ^ j>e...f 

c "' d dx e dxf dx c dx d "' k 
A tensor can be defined at a single event, or along a curve, or on the whole of space-time, 
in which case the components are functions of the coordinates and we call T a tensor field. 
If q = then T ia a contravariant tensor; if p = 0, it is a covariant tensor. 
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An object that behaves as a tensor under change of local inertial coordinates at an event 
determines a tensor at the event under general coordinate transformations. 
Note that since 

d¥dx~ b=6 ^ (24) 
one could equally well write the transformation law with all the tilde ( ~ ) and un-tilde 
quantities interchanged. 

Example. The metric g a b is a tensor field of type (0, 2). It has the transformation law 

9ah= ~ 9cd dx~^- (25) 

Example. The contravariant metric has components g ab , where (g ab ) is the inverse matrix 
to g a b (i.e. g ab gbc = It is a tensor of type (2,0). This is proved from (25) by the 
following steps, which are well worth following carefully because they illustrate some basic 
techniques of index manipulation. The proof makes several uses of (24). 

Qx b 

(multiply both sides by — — and sum over b) 

ox e 

(multiply by g a ^g eh and sum over a, e) 
(multiply by — - and sum over f ) 

OX J 

Example. The gradient d a f of a scalar function is a covector field (tensor of type (0, 1)). 

Example. If x a = x a (s) is the worldline of a particle parametrized by s, then V a = dx a /ds 
is a 4-vector field along the worldline. If g a bV a V b = g a \,x a x b = 1, then V is called the 4~ 
velocity and s is called the proper time. 

We can carry out all the operations on tensors in exactly the same way as in SR 
with the exception of differentiation: the components d a T d... no longer obey the tensor 
transformation law. We shall come back to that. 

Indices are raised and lowered by contracting with g ab and g a b, but this now involves 
more than just changing the signs of a few components. For example, if T ab c is a tensor of 
type (2, 1), then the contraction T ab b is a tensor of type (1, 0) (one free upper index a). If 
ot a is a covector then g ab ot c is a a tensor of type (2, 1) and its contraction 

a a = g ab a b 



dx b 


dx c 


9ab dx* 


9ce dx a 


9x1 eh 

- — g en 
dx e 


_ ~ g af dxh 

dx a 


g kh 


„ f dx h dx k 
dx a dxf 
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is a vector. This is the operation of raising the index. One similarly lowers indices, for 
example by putting X a = g a bX b . Raising followed by lowering returns to the starting point 
since 

g ab g bc = 8 c a . 

The exceptional operation, differentiation, is more subtle in a general space-time. In the 
next lecture, we start by considering how to define a derivative operation for vectors. 



6.4 Summary of the mathematical formulation 



Space-time is a four-dimensional manifold M with a metric tensor g a t, (symmetric, type 

(0,2), signature H ). 

The points of M are the events. If x a = x a (u) is the worldline of a particle, where u is 
a parameter, then 

dx a dx b 

g a b-;—-j—<iu 

du du 

is the proper time along the worldline: the time measured by a clock carried by the particle 
(this the clock hypothesis) . The 4- vector V with components V a = dx a /ds is the 4- velocity. 
The metric determines the behaviour of free particles via the geodesic hypothesis 



dV 



ds 2 



+ r? 



6c 



dx a dx b 
ds ds 







where s is proper time. 

If A is an event, then there exists a local coordinate system such that x a 



at A and 



(9ab(x)) 



(1 








o \ 





-1 














-1 





Vo 








-v 



0(2). 



In these coordinates, T^ c = at the origin (the event A). Such a coordinate system is 
interpreted as the local inertial coordinate system set up by an observer in free fall at A. 
We identify 4-vectors and tensors at A with vectors and tensors in SR by taking their 
components in local inertial coordinates. 

The metric determines an inner product g(X, Y) = X a Y a on the space of 4-vectors at 

an event with signature H (that is, it is not positive definite). As in SR, we say that 

X is timelike if X a X a > 0, and so on. 
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Lecture 7 



7.1 Differentiation of vectors 

We need to be able to differentiate vectors and tensors — for example to be able to write 
down Einstein's equation's. But there is a problem — illustrated by the following examples. 

Example. Let X a be a vector field and put d b = d/dx b . Then 

*** = sS/nOn* 

ox ox c \ ox a 
dx c dx a q dx c d 2 x a 
dx b dx d c dx b dx c dx d 

The first term is what we want for a tensor transformation law. The second term is the 
problem (in SR, where the coordinate transformations are all linear, it vanishes automati- 
cally: the problem is that we must now consider general, nonlinear coordinate transforma- 
tions). 

Example. Let x a = x a (s) be the worldline of a particle parametrized by proper time. 
Then 

dx a dx a dx b 



ds dx b ds 

which implies that the 4-velocity components transform in the right way. But 

d 2 x a dx a d 2 x b d 2 x a dx b dx c 
+ 



ds 2 dx b ds 2 dx b dx c ds ds 

Again the second term is the obstruction to a nice transformation law: the obvious defini- 
tion of 4-acceleration does not lead to a tensor. 

There is however a way out that does lead to a tensor transformation law in both these 
cases. Recall that the Christoffel symbols 

r fe c = y ad (d b gdc + d c g ba - d a g bc ) 

obey the transformation law 

dx a dx e dx? ~ , dx a d 2 x d dx a dx e dx? ~ , d 2 x a dx e x? 

j^a T^a t^u 

1 hr — «~\ — j ~c\ r a r 1 P f t ~_ j 7: r^. : — ^~ j 7; :t; z l P f ~ 



bc dx d dx b dx c e f dx d dx b dx c dx d dx b dx c e ^ dx e dxf dx b dx c 
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The second term is exactly what we want to cancel the unwanted term in the first example. 
Therefore 

V b X a = 3 b X a + T a bc X c 

transforms as a tensor of type (1,1). It is called the covariant derivative of X a . 
To show the transformation rule in detail, 



dx c d f dx h ~„ A dx b dx$ dx h ^ , rA d 2 x b dx e dx? 
dx c dx b 



a cd ~ dx a dx c \dxf dx e dx c dx d fh dx e dxf dx c dx d 



d c X d + r a ,x< 



dx a dx d \ c de ' 

In a coordinate system such that d a g bc = at the event x a = 0, we have T^ c = 
at x a = and hence V a X b = d a X b . Note that this holds only at the origin. We could 
have used this property to define the covariant derivative: that is, we could equally well 
have defined the covariant derivative by requiring that the value of V a X b at A should be 
the tensor that coincides with d a X b in local inertial coordinates at A. Then the tensor 
transformation law would have enabled us to write down its components in a general 
coordinate system. It is a common technique to define a tensor by giving its components 
in a particular coordinate system and then to use the transformation law backwards. 



7.2 Parallel transport 

Consider two nearby events A and B with coordinates x b and x b + 5x b . To the first order 
in 5x b , 

Sx b V b Y a = 5x b d b Y a + 5x b Y a bc Y c = Y a (x + 5x)-(Y a (x) - 5x b T a bc Y c (x)^j . 

Thus the covariant derivative compares Y a (x + 5x) (the components of Y at B) with 
Y a (x) — 5x b T bc Y c , which we think of as the result of displacing Y a to the 'most nearly 
parallel vector at B\ 

Definition 7 

The vector at B with components Y a (A) — 5x c T bc Y c (A) is said to be obtained by parallel 
transport of Y a from A. 

In local inertial coordinates at A, F = at A, the vector at B is the one with the same 
components (to the first order in Sx). It makes more sense to express this in terms of 
parallel transport along a curve — we then don't have to worry about infinitesimals. 



34 



Figure 8: Parallel transport of Y along the curve x a = x a (u). 

Definition 8 

A vector Y is a parallel transported (or propagated) along a curve x a = x a (u) whenever 

du +lbc du " U - 

This is a set of ordinary differential equations for the components Y a as functions of the 
parameter u: it determines the Y a in terms of their values at the initial point of the curve. 
The tangent vector to a geodesic, which is a solution of 

x a + T a bc x b x c = 0, 

is parallel propagated along the geodesic: this is the sense in which geodesies are curves in 
a curved space-time which are 'as straight as possible'. 

Parallel propagation around a closed curve need not bring you back to your starting 
point: this is a manifestation of curvature. 



7.3 The wave equation 

Suppose that u is a function on space-time, 
ponents of a covector, the gradient covector. 
V a M by putting 

V a w = 



Then the partial derivatives d a u are the com- 
We can define a vector field with components 

g a %u. 
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That is, by raising the index. This is the gradient vector. The wave operator sends u to 

V a (V a u) = d a (g ab d b u) + T a ab g bc d c u . 
It is an exercise on problem sheet 4 to show that 

1^=100 log 

where g is the determinant of the matrix (g a b) — it follows from the fact that if A is a square 
matrix depending on the coordinates x a , then 

<9 a logdetA = tr(A~ 1 d a A) . 

This implies that 

r a& = y M (db9ad + d a g bd - d d g ab ) = \g bd d a g bd = d a log ^/\g~\ . 
Thus the wave operator is also 




It is a natural operator; that is, it is independent of the choice of coordinates. 
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Lecture 8 



8.1 Covariant derivatives of tensors 

We extend the defintion of V a to scalars by putting V a / = d a and to covectors by 

V a a b := d a a b - a c T c ab . 

By a similar argument to that used in the case of vectors, this transforms as a tensor of 
type (0,2). 

Exercise. Show that d a (a b Y b ) = (V a a b )Y b + a b V a Y b . (Note that a b Y b is a scalar, so the 
gradient covector on the LHS is well defined.) 

For a general tensor field, the covariant derivative has one gamma for each index. For 
example, 

W a T bc d = d a T bc d + T b ae T ec d + K e T be d ~ K d T bc e . 

We also write V a / for d a f when / is a scalar. 
Properties of the covariant derivative. 

(cdl) V a (T~... + S-J = V„7- .. + V a S-.... 

(cd2) V a (/T-...) = ./T„7 + (V a /)T-.... 

(cd3) V a (T-..S-J = V a (T-„.)S-... + T~.y a (S~J. 

(cd4) V a T bc b is the same whether the contraction is done before or after the 
differentiation. 

(cd5) The covariant derivative of the Kronecker delta vanishes since 

VJ b c = d a 5 b + T b ad 6 d c - T d ac 6 b d = 0. 

(cd6) For a scalar /, but not for a general tensor, 

v a v 6 / = d a d b f - djr c ab = v b v a f. 

(cd7) The covariant derivative of the metric tensor vanishes, since 

V n fc = d a g bc — gdc^ih ~ 9bdF d ac 

= d a g bc - \{d a g cb + d b g ac - d c g ab } - \{d a g bc + d c g ab - d b g ac } 
= 
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(cd8) V a g bc = 0. This follows from (cd.7) and 

= V«($) = V a (g bc g cd ) = V a (g bc )g cd + g bc V a g cd = 0. 
Hence V ' a g bc = 0. 

It follows from (cd7) and (cd8) that raising and lowering can be interchanged with covariant 
differentiation. For example if X a is a vector field, then V a X b is well-defined: it does not 
matter whether you lower the index on the X before or after the differentiation. 

Example. Maxwell's equations. In a curved space-time and in the absence of sources, 
these are 

V a F ab = 0, V a F bc + V b F ca + V c F ab = 

because these equations are covariant and reduce to the special relativity form in local 
inertial coordinates at a point. Gravity affects light through the Ts. 

8.2 Connections 

All that is needed to define the covariant derivative in a coordinate independent way is the 
transformation rule 

dx a ~ j dx e dx? dx a d 2 x d 
bc ~ dx d ef dx b dx c dx d dx b dx c 
A field of r^ c s with this transformation property is called a set of connection coefficients. 
The corresponding operator V is called a connection — it connects the spaces of vectors and 
tensors at nearby events. 

(cdl)-(cd5) hold for all connections; (cd6) holds only if the connection is torsion-free — 
that is Y c ab = T ba . (cd7) holds in addition only for 

Kc = h ad (d b g cd + d c g bd - d d g bc ). (26) 

This is the unique torsion-free connection for which the covariant derivative of the metric 
tensor vanishes (problem sheet). It is called the Levi-Civita connection. 

Note that if T bc is a set of connection coefficients and Ql c is a tensor, then r^ c + Ql c is 
also a set of connection coefficients. All connections can be obtained in this way once once 
is given. 

From now on V will always be the Levi-Civita connection, defined by (26). 
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8.3 Curvature 



In Minkowski space, there are global coordinate systems in which g ab is constant: in such 
coordinates V a = d a and therefore V a V& = VbV a (acting on vectors). So if in a general 
space-time, V a 7^ V b V a on vectors, then we know that the metric cannot be reduced to 
the SR form by a coordinate change. 

Proposition 6 

For any metric g ab , there is a tensor field R abc d of type (1, 3) such that 



V a V b X d - V b V a X d = R abc d X c 



for any 4- vector field X. 



The tensor R ab d is called the Riemann tensor or the curvature tensor. 



Proof. From the definition of V, 

V a V b X d = V a (d b X d + T d bc X c ) 

= d a d b X d + (d a T d bc )X c + T d bc d a X c 

+ T d ae (d b X e + Tt c X c ) - Tl b {d e X d + T d ec X c ) 

Hence 

(V a V fe - V b V a )X d = (d a T d bc - d b T d ac - T d be T e ac + T d ae r bc )X c . 

We define the expression in brackets to be R ab d - We must show that it is a tensor. The 
direct method is horrible. We know, however, that the left-hand side is a tensor. Hence, if 
we change coordinates, 

~ ~ ~ ~ ~ ~ Q^xP dcc^ doc^ 

v a v b x d -v b v a x d = Wa Q^ b Q^R m c s x c 

dx p dx q dx r dx d s ~ c , . 

W a lB b dx~ c ~dx~ s pqr ' ^ ' 

However, had we worked from the beginning in the tilde coordinates, we would have ob- 
tained 

(V Q V 6 - V b V a )X d = R ab d X c , (28) 

where R ab d is defined in the same way as R ab d . Since (27) and (28) hold for any X, we 
deduce that 

~ d dx p dx q dx r dx d s 

abc = W a ^x~ b lB c dx~ s pqr ' 
which is the tensor transformation law. ■ 
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Corollary 7 

If there exists a vector field X such that V a V&X d ^ VbV a X d , there does not exist a 
coordinate system in which the metric coefficients are constant. 



8.4 Symmetries of the Riemann tensor 

In terms of the connection coefficients 

p d a -pd o -pd -pd -pe , -pd -pe 

K-abc — u a L be u b L ac 1 be L ac 1 ae L be 

Pick an event A and choose coordinates such that d a gbc — at A. Then T^ c = at A. We 
can drop all the first derivatives of g a b and g ab , but not the second derivatives. So, at the 
event A (but not elsewhere), 

Rabcd = gdedaiXbc) ~ 9dedb(Tac) 

= \d a {d c gbd + d b g dc - ddgbc) - \d b {d c g da + d a g dc - d d g ac ) 
= \ [d a d c gbd + dbddg ac - d a d d g bc - d b d c g ad ] 

From this we deduce that 



(SI) 


Rabcd 


Rbacd 


(S2) 


Rabcd 


Rcdab 


(S3) 


Rabcd 


Rabdc 


(SA) 


Rabcd 


+ Rbcad + Rcabd = 



Note that symmetries (SI), (S3) and (S4) imply (S2) (exercise). 
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Lecture 9 



9.1 Bracket notation 

For a general covariant tensor with p lower indices, 

T[ab...c] ■= -j si g n (° r ) T ( 7(a) CT (b)... CT ( C ) 

perms 

T(ab...c) ■= -j- T ^Mb)...<x( C ) 

^' perms 

where the sums are over the permutations a of p objects, and sign (a) is 1 or —1 as a is 
even or odd. For example, 

T[ab] = \{T a b — Tb a ) 

T(ab) = \{T a b + T ba ) 

T^abc] g (T a bc T^bca T^cab T bac T^acb -^cfea) 

-^(afec) g (-^afec "I - -^~ftca "I - Tcab Tbac -^acb -^~c6a) 

The same definitions apply to brackets on a subset of the indices and to brackets on upper 
indices. For example 



rp[ab](cd) l^rpabcd rpbacd _|_ rpabdc rpbadt , 



There is a possibility of ambiguity over the order of the operations if two sets of brackets 
overlap (as, for example, in the expression T[ (& c ]d)); so overlaps are forbidden. Nested 
brackets, however, are unambiguous, although they can always be simplified since 

7 [-(-)...] = = 7 (...[..:...) ' = 7 t-] ' 7 (-(...)...) = T (-..) • 

Examples. With this notation, the fourth symmetry of the Riemann tensor reads 

R[abc]d = . 

Maxwell's equations (without sources) are 

V a F ab = 0, V [a F bc] = . 

The symmetries of the contravariant metric g ab and of the alternating tensor e a b c d (see 
problem sheet 1) can be expressed respectively as 

g[ab] _ tghd = e[abcd] ■ 
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9.2 The Bianchi identity 

Choose coordinates such that T bc = at an event. Then 

V a R bcd e = d a d b T e cd - d a d c T e bd + terms in TOT and TIT 
Since the first term on the RHS is symmetric in ab and the second in ac, 

at the event in this coordinate system. However, this is a tensor equation, so it is valid in 
every coordinate system. It is called the Bianchi identity. 

9.3 The operator D 

Let w be a the worldline of an observer (not necessarily in free fall). Let s be the proper 
time along uj and let V a = dx a /ds be the 4- velocity of O: this is a vector defined at each 
point of uj. 

We need a tool: a derivative operator that measures the rate of change of vectors and 
tensors as observed by O. Suppose that Y a is a vector field. Then at each event on uj, we 
define DY b by 

dr a dY b 
DY b = V a V a Y b = —d a Y b + T b V a Y c = — + T b V a Y c . 
ds ds 

The first equality makes it clear that DY a is a well-defined vector at each point of u. The 
last, that the values of DY a along uj depend only on the values of Y a (s) along u: so DY a 
makes sense for vectors fields that are defined only along uj. Note that DY a = is the 
equation of parallel transport. 

We similarly define D on tensor fields along uj by, for example, 

dT a 

dt: = ^ + r a cd v c T d b - T d cb V c T a d . 

If the observer is in free fall, then D = d/ds in local inertial coordinates in which the 
observer is at rest at an event. 
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9.4 Geodesic deviation 

The crucial property of curvature is that it measures the relative acceleration of nearby 
particles in free fall. That is, it directly encodes the gravitational field. 

Imagine a cloud of particles in free fall. Let us suppose that an observer O is on one 
of the particles, and that this particle has worldline u>. Suppose that the observer looks 
at a nearbly particle and measures its position in local inertial coordinates. In SR, it will 
move in a straight line at constant speed, and will have no acceleration. What happens in 
a gravitational field? 

The 4- velocities of the particles form a vector field V a . Since the individual particle 
worldlines are geodesic, 

dV a 

V a V a V b = DV a = — + TlV b V c = 0. 
ds 

Pick out a particle P near O, and at each event on u, let Y a be the 4- vector joining 
the event to a simultaneous event at P (since P is 'near' O, Y is small: we shall ignore 
second order terms in its components). In the LIC in which O is instantaneously at rest, 

Y has components {0,y), where y is the position of P. If u is given by x a = x a (s) in 
general coordinates, P's worldline is x a (s) + Y a (s) + 0(2), where 0(2) denotes second- 
order and smaller terms in the components of y. The parameter s here is the proper time 
along the worldline of O. The proper time separation ds between two nearby events x a (s) 
and x a (s + 5s) on the worldline of O is the same to the second order in y as the proper 
time between the corresponding events, x a (s) + Y a (s) and x a (s + 5s) + Y a (s + 5s), on 
the worldline of P. To this approximation, therefore, s is also the proper time along P's 
worldline. 

We note that Y a is a vector field along u such that V a Y a = (because Y = (0, y) and 

V = (1, 0) in the LIC in which O is instantaneously at rest). 

Since DV a = 0, we have 

= D(V a Y a ) = V a DY a and = D(V a DY a ) = V a D 2 Y a . 

In the local rest frame of O at an event on its worldline, the 4- velocity of O is (1,0); and 
the vectors Y a , DY a and D 2 Y a are, respectively, of the form (0, y), (0, u) and (0, a), where 
u is the relative velocity of P to O and a is the relative acceleration. 

We are interested in the relative acceleration, and therefore in D 2 Y a . We want to 
express this in terms of the curvature. The key to this is the following result. 
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Proposition 8 
DY a = Y b V b V a . 

Proof. We know that 



V a (P) = — + — + 0(2) = V a (0) + -^ + 0(2) . 
ds ds ds 

On the other hand, by expanding to the first order in e, 

V a {P) = V a {0) + Y c d c V a + 0(2) . 

Therefore dy a /ds = Y c d c V a . It follows that 

dY a 

jjya = + r aybyc = ycQ ya + pa ybyc = y 6y ya 

ds 

which is the result we need. ■ 

Now we can derive the equation of geodesic deviation (Jacobi equation), which is of 
central importance to the physical interpretation of curvature. 

D 2 Y d = D(Y b V b V d ) 

= (DY b )V b V d + Y b D(V b V d ) 

= (Y a V a V b )V b V d + Y b V a V a V b V d 

= Y a (V a V b )V b V d + Y b V a V b V a V d + R abc d V a Y b V c . 

But 

v a v b v a v d = v b (v a v a v d ) - (v b v a )(v a v d ) = -iy b v a )iy a v d ) 

since V a, VaV d = by the geodesic equation. Therefore the first two terms in the last line 
cancel, and 

D 2 Y d = R abc d V a Y b V c , 

which is the geodesic deviation equation. It relates the relative acceleration of nearby 
particles in free fall to the curvature tensor. 
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Lecture 10 
10.1 Tidal forces 

The relative acceleration of two particles in free fall is determined by the equation of 
geodesic deviation 

D 2 Y d = R abc d V a V c Y b . 

If you are travelling with one of the particles, with 4-velocity V, then the acceleration a of 
your neighbours is a linear function of their position. In local inertial coordinates in which 
you are (instantaneously) at rest, V = (1,0), Y = (0,y), where y is the position vector of 
a nearby particle, and 

a = My 

where M is the 3x3 matrix with entries 

M ij = ROW 

= 1,2,3). It is symmetric because of the symmteries of the Riemann tensor, and 
because the compoents are taken in local intertial coordinates. 

What is the corresponding result in Newtonian gravity? Consider a cloud of particles 
in free fall. The acceleration of each particle is given by r = — V0, or in components, 
by r\ = —di(p. Since the particle at O has acceleration (—d^)o and the particle at P has 
acceleration (— di<f))p, the acceleration of that at P relative to that at O has components 

Oi = {-d i( j))p - {-d4) = -yAdit + 0(2) (29) 

to the first order in y, where y is the position vector of P relative to O, the second 
derivatives are evaluated at O, and there is a sum over j = 1, 2, 3. 

Aside. Another way to read (29) is in terms of tidal forces. If O and P are joined by 
a light rod, and falling towards the moon with the rod pointing towards the moon, then 
then there is a tension force in the rod. This is true whichever particle is in the lead. If we 
think of O as the centre of the earth, and of P as a mass of water on the surface, then we 
see that this tidal force will push P away from the centre: this is true whether P is on the 
surface directly under the moon or on the opposite side of the earth. Thus the tidal force 
raises two 'humps' in the ocean, one under the moon and one on the other side. As the 
earth rotates, the humps move round, giving two high tides each day. We can trace the 
reason that there are two high tides to the symmetry of the second derivatives in djdi<f). 
Before Newton, even Galileo's explanation of this fact was confused (and erroneous). 
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When our observer in curved space-time looks at the nearby particles, therefore, he 
reckons that he is in a gravitational field with potential 4> such that 

Mij = -didjcj) = R Qi0 ° 

at his location. Now in empty space, Poisson's equation reduces to V 2 = 0. That is 
didi4> = or, in other words, trM = 0. So in GR, we require R ab b V a V c = 0. Since this 
must hold for every 4-velocity V, we are led to Einstein's vacuum field equation 

Rab = 0, (30) 

where R a b is the Ricci tensor, defined by R ab = R acb c . The vacuum equation is in fact 
ten equations (for the ten independent components of the symmetric tensor R ab ) in ten 
unknowns (the ten independent components of the metric g a i). The equations are nonlin- 
ear, as anticipated. The use of the summation convention makes them look very simple. 
Written out explicitly without this notation, each equation would contain 1248 terms. Not 
surprsingly, therefore, it is not easy to find solutions. 

We shall see two justifications for the vacuum equation: (1) it reduces to the Newtonian 
equation in the weak field limit, and (2) it has a solution (the Schwarzschild solution) like 
<p = —GM/r, which encodes the inverse-square law of gravity. 

10.2 The weak field limit 

We consider the case in which the gravitational field is weak and slowly varying. That is 
the metric is 

Wab = gab + th ah , 

where g a b = diag(l, —1, —1, —1) is the Minkowski space metric, e is a small parameter, and 
the coordinates are the Minkowski space inertial coordinates t, x, y, z. Weak' means that 
we ignore terms of order e 2 ; 'slowly varying' means that doh a b = 0(e). 

To find the vacuum equation in this case, we first have to find the contravariant metric 
w ab (defined by w ab Wb c = 5"). This is given by 

w ab = g ab + O(e) , 

where g ab = diag(l, — 1, — 1, — 1) is the contravariant Minkowski space metric . There is 
an immediate possibility for confusion here since we are dealing simultaneously with two 
metrics, w and g, so we shall avoid raising and lowering indices. 
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The Christoffel symbols are given to the first order in e by 

T a bc = \w ad {d b w cd + d c w bd - d d w bc ) = \eg ad {d b h cd + d c h bd - d d h bc ) + 0(e 2 ) 

and therefore the Riemann tensor is 

R-abc — 9 a T d c - d b T d ac + 0(e 2 ) 

= \eg de (d a d c h he + d b d e h ac - d a d e h bc - d b d c h ae ) + 0(e 2 ) . (31) 

Now consider the motion of a slow moving particle with worldline x a = x a (t). Since 
the speed u is small, we have 7(w) ~ 1 and so we can identify the coordinate time t with 
the proper time along the wordline. The 4-velocity is 

(V a ) = (^j = (1, 0, 0, 0) + 0(6), s = t + 0(e) , 
d 2 x a ^ n dx b dx c 



and the geodesic equation is 



To the first order, this is 



where 



ds 2 bc ds ds 
d 2 x a ^ n 

J- nn 



ds 2 ■ -°° 



0(e 2 ) when a = 



r °° \\td a h m + 0(t 2 ) when a ^0, 

since doh ab = O(e). When a = 0, this gives us d 2 t/ds 2 = 0(e 2 ). The other three compo- 
nents give the equation of motion 

r = -lV(eh 00 ) + O(e 2 ), 

where the dot can be differentiation with respect to either t or s. Thus if we want to reduce 
GR to the Newtonian theory in this limit, then we must take <ft — ^h 00 . 

The derivatives of h ab with respect to x° in (31) are all of order e. So we have 

Row" = leg de d b d e h o + 0(e 2 ) . 

Therefore the 0, component of the vacuum equation R ab = is 

Row — l e (-9 2 h 00 - d%h 00 - dlh 00 ) = -\eV 2 h 00 = , 

that is, V 2 = 0, which is the Laplace's equation (the vacuum equation in Newtonian 
theory). Thus in this limit, Einstein's theory reduces to Newton's. (What about the other 
nine components of the vacuum equation?). 



47 



10.3 The non- vacuum case 

What happens when there is matter present? What is the analogue of Poisson's equation 
V 2 = 47rGp? 

We shall consider only the case in which the matter generating the gravitational field 
is a dust cloud. Then its energy density is encoded in the energy-momentum tensor T ab = 
pjjajjb w h ere jja j s 4-velocity field of the dust and p is the energy (mass) density 
measured in the local rest frame. We know that d a T ab = in local inertial coordinates since 
this continuity equation holds in SR. Therefore we have V ' a T ab = in general coordinates, 
since this is a tensor equation which reduces to the continuity equation in local inertial 
coordinates. This suggests that the GR field equation should equate R ab to a constant 
multiple of T ab . Unfortunately, this will not do because in general V a R ab 7^ 0. But 
there is a tensor closely related to the Ricci tensor which can be put on the LHS without 
contradiction: this is the Einstein tensor 

G a b — Rab — \RQab , 

where R = R a a is the Ricci scalar or scalar curvature. 
Proposition 9 

For any space-time metric, V a G ab = 0. 
Proof. The Bianchi identity is 

V ' a Rbcde + ^bRcade + V ' c R a bde = , . 

By contracting with g ad g ce , we obtain 

= 2V a R ab - V b R = 2V a (R ab - \g ab R) = 2V a G ab , 

which completes the proof. ■ 

Our candidate for the field equation is G ab = \pU a U b , with A constant. By contracting 
with g ab , we obtain 

R-2R = AT/ = Xp 
since g ab g a b = 4 and U a U a = 1. So an equivalent form of the equation is 

R a b = A(T afe — \pg a b)- 
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Now in the coordinates we used in the weak field limit, R 00 = — V 2 0, while T 00 = p. 
Thus in this limit, we have V 2 </> = — ^Ap. To obtain the correct correspondence with the 
Newtonian theory, therefore, we must take A = — 8nG, which means that the field equation 
is 

Rab — 2^-9 ab = ~ 87T GT a b . 

This is in fact quite general: the tensor T a b can be the energy-momentum tensor of any 
form of matter or energy. 

In future, we shall always use units in which G = 1 = c (note that given the unit of 
time (say the second), putting putting c = 1 fixes the unit of distance (the light-second) 
and putting G = 1 fixes the unit of mass. 

Example. Calculate the conversion factors to SI units and take note that the units we 
shall use are not likely to find favour with engineers. 
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Lecture 11 



11.1 Spherical symmetry 

The next few sections will be of the form an extended worked example. 

We want to find the gravitational field outside a spherical body of mass m, that is the 
solution of the vacuum equation analogous to the Newtonian potential <fi = —Gm/r (for 
simplicity, we shall now choose units so that G — 1). By saying the solution has mass m, 
we mean that goo ~ 1 — 2m /r for large r, so that a long way from the body, where the 
field is weak, it looks like the field of a static spherically symmetric body of mass m: in 
operational terms, m is the mass measured by analysing orbits in the field of the body 
near infinity. 

We want the metric to have the symmetries appropriate to a static spherical body. Now 
in spherical polar coordinates, the Minkowski space metric is 

di 2 - dr 2 - r 2 (d6 2 + sin 2 6d<f) 2 ) . (32) 

The expression in brackets is the metric (first fundamental form) on the unit sphere. 

Our metric must reduce to (32) when m = and in any case in the limit r — > oo. This 
flat metric has the following features. 

• No t-dependence. 

• No dtdr, dtd0, or dtdd terms: it is therefore time reversible — that is, invariant under 
t i-> -t. 

• No drdd or drd(p terms: thus at constant time, the radial vector is perpendicular to 
the surfaces of constant r. 

• The metric on the surface t = const, r = const is a constant multiple of the metric 
on the unit sphere. 

• The coefficients of dt 2 and dr 2 are independent of 9 and 4>. 

The last three properties might be taken to characterize spherical symmetry. We shall 
assume that out metric has all these properties, and thus that it is of the form 

ds 2 = A{r)dt 2 - B{r)dr 2 - C(r)r 2 (d6 2 + sin 2 6d<p 2 ) , 

for some functions A, B, C of r. We can in fact take C — 1 without loss of generality by 
replacing r by r\/~C ' . Our task then is to solve the vacuum equation R a b = 0, subject to 
the boundary conditions A, B — > 1 and A = 1 — 2m /r + 0(r~ 2 ) as r — > oo. 
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11.2 The curvature tensor 

We need to find R a b in terms of A and B (with C — 1). We shall do this by a straightfor- 
ward, but head-on method: a more elegant method can be found in Hughston and Tod, 
but it requires some understanding of differential forms. 

The first step is to find the Christoffel symbols from the geodesic equations. These are 
the Lagrange equations 



with x° = t, x l = r, x 2 = 8, x 3 = (p. The idea is to read off the Christoffel symbols by 
comparing the Lagrange equations with x a + T^ c x b x c = 0. 
The Lagrange equations are 




generated by the Lagrangian 



L = \{Ai 2 - Br 2 - r 2 6 2 - r 2 sin 2 6<p 2 ) , 




















where the dot denotes the derivative with respect to s. Equivalently 



f + 



lA'B-H 2 + \B'B- X r 2 



i+A'A'Hf 
- B~ 1 r9 2 - B^r sin 2 6 <p 2 
9 + 2r _1 0r -sin^ cos eft 
<p + 2r~ 1 (j)f + 2cot 9 6(f) 















0. 



From these we read off the Christoffel symbols 
b 7^ c: where do they come from?) 



r^ c (note carefully the factors of | when 



(o = 0) 




1 ii 



A'/2A 
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(a = l) Fl = A'/2B, Fl^B'fiB, Y\ 2 = -r/B, Y\ z = —r sin 2 9/ B 
(a = 2) r 2 .! = r 2 2 = r" 1 , r 2 3 = -sin^cos^ 
(a = 3) rj^ = T? 3 = r" 1 , 113 = 1*= cot0. 

All the others vanish. From the definition of the curvature tensor, we have 

p d Q -pd Q p-d -pd pe , -pd pe 

-n-aftc — ^a 1 be °b L ac 1 fee 1 ac 1 ae 1 fee ■ 

The components -R ac of the Ricci tensor are then given by putting b = d and summing. 
We find 

-^232 3 = ^32 - d 3 T\ 2 - Y^T^ + r^r^ 

= d e (cot9) + B- 1 + cot 2 9 

= -l + B- 1 
R^ 1 = -A"/2B + B'A'/AB 2 + A 2 /ABA 
R 2o = Row = ~A'/2Br 
Rni = Rni = ~B'/2Br 
R W i = —BR 01 q I A 

R 3Q3 ° = -r 2 sin 2 9R J /A = r sin 2 9A'/2B A 
R^ 1 = r 2 sin 2 9R 13l 3 /B = -r sin 2 9B'/2B 2 . 

Finally, we obtain the vacuum equations in the form 

Roo = R 01 o + R 2o + R 030 = ~A"/2B + B'A'/AB 2 + A' 2 /ABA - A' /Br = (33) 

together with 

R n = A"/2A - A' 2 /AA 2 - B'A'/ABA - B'/Br = (34) 
#22 = R 33 / sin 2 9 = rA'/2BA - rB'/2B 2 + 1/5-1 = 0. (35) 

All the other components of the Ricci tensor vanish identically (as can be seen by direct 
calculation or by using the fact that R a t, must have the same symmetries as the metric). 

In all, we have three equations in the two unknowns A, B. Fortunately they are 
consistent. 
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If we take A times (33) and add A times (34), we get 

AB' + BA' = , 

and hence that AB is constant. Since we want A, B —>■ 1 as r —>■ oo, we must therefore 
have AB = 1. By substituting into (35), we then get that rA' + A = 1 and hence that 

1 k 
A = - = l + - 
B r 

for some constant k. But for large r, we want A = 1 — 2m/r + 0(r 2 ), so k — —2m, and 
the solution is 

dt?=(l-—) dt 2 d ^—- r W + sm 2 9d ( f> 2 ). 

\ r J 1 — 2m/r 

This is the Schwarzs child metric. The method of derivation is notable only for the incentive 
it gives to find more subtle methods for tackling Einstein's equations. 
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Lecture 12 



12.1 Stationary observers 

An observer in a fixed location (relative to our coordinate system) has a wordline with 
constant r, 9, <p, and therefore has 4-velocity U with only the first component nonzero. 
Since U a U a = 1 and U° > 0, therefore, we must have 

(U a ) = (1/0 -2m/r, 0,0,0). 

This worldline is not geodesic — as we know, for example, from the fact that an observer at 
rest on the earth's surface is accelerating relative to the local inertial frame. He interprets 
this acceleration as the 'force of gravity'. 

In local inertial coordinates at an event on his worldline, the four- acceleration of the ob- 
server is a a = dU a /ds; in arbitrary coordiantes, therefore, a a = U h VbU a . The acceleration 
actually felt by the observer is y/—a a a a 

By using the calculation of the Christoffel symbols in the last lecture, we have 

(U b V h U a ) ={U°d U a + T a 00 U°U°)= (0, A'{U°) 2 /2B, 0, 0) = (0, \A' , 0, 0) 

where A = 1/B = \—2mjr). Thus the 4-acceleration of the observer is (a a ) = (0, m/r 2 , 0, 0), 
as one might expect by naive analogy with Newtonian theory. However, the acceleration 
felt by the observer is 

g = V~a a a a = ——==== (36) 
ryl — 2m jr 

Thus the 'force of gravity' is given by the same inverse-square law g = m/r 2 as in Newtonian 
theory for large r, but increases to infinity as r approaches the Schwarzschild radius r = 2m. 
We shall see later that r = 2m is the 'event horizon' of a black hole. What we are observing 
is a consequence of the fact that inside a black hole, you would have to travel faster than 
light in order to stay in the same place (what would the Red Queen say?). 

12.2 Potential energy 

The t equation for the geodesic motion of a free particle is 

d fdL\ 
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since the Lagrangian is independent of t. Therefore E — (1 — 2m/r)t is constant along the 
particle worldline. What is the interpretation of this constant? 

Suppose that the particle has 4-velocity V and unit mass. Then relative to our observer 
'at rest' at some point in the particle's history, the particle has speed v given by 

7(u) = -=L= = U a V a = g 00 U°V = iy/l-2m/r. 
V 1 — v z 

Therefore 

a/1 - 2m/r 

Cj = . . 

For large r and small v, this is approximately 

E = 1 + |u 2 — m/r + smaller terms . 

Thus is the sum of the rest energy (Mc 2 with M = 1 and c = 1), the kinetic energy |t> 2 
relative to the observer, and the Newtonian potential energy —m/r. Thus it is reasonable 
to interpret E as the total energy of the particle. We note that this is consistent with (36), 
which can be written g = d r *Jl — 2m/ r, and which therefore also suggests that we should 
interpret v/l — 2m /r as the potential energy of a unit mass particle at rest. Conservation 
of energy is a consequence of dL/dt = 0. 

12.3 Photons 

Geodesies in space-time are generated by the Lagrangian 

L = ^9abX a £ b ■ 

In the case of free particles, the dot denotes differentiation with respect to proper time 
s — the time measured by a clock carried by the particle. In the case of photons, the dot is 
differentiation with respect to an affine parameter, which is defined only up to a constant 
factor (and of course addition of a constant). 

In Minkowski space, a photon worldline is a null line. The frequency 4-vector K is 
tangent to the worldline and encodes information about the frequency of the photon, 
as measured by a moving observer. If the observer has 4-velocity U, then the observed 
frequency is uj = U a K a . 

In special relativity, the frequency 4-vector is constant along the photon worldline. By 
our usual principle that special relativity should hold over short times and distances in local 
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inertial coordinates, in general relativity, we must have that K is tangent to the photon 
worldline (now a null geodesic) and that K a V a K b = 0. But if we put W a = dx a /da, then 
the geodesic equation is W a 'V a W b = 0. So W is a constant multiple of K, and by rescaling 
a, we can take W = K. Then the frequency 4- vector is given by K a = dx a /da. 



12 A Gravitational redshift 

In the Schwarzschild metric, 



r J 1 — 2m I r 



(37) 



and (1 — 2m/r)t is constant along both timelike and null geodesies. 

Consider two observers 0\ and 2 at rest relative to the Schwarzschild coordinates at 
r = r 1 and r = r 2 . If 0\ sends out a photon to 2 , and if the frequency measured by Oi 
at transmission is Ui, then what is the frequency at reception as measured by 2 ? 

Let the photon's worldline be x a = x a (a), where the affine parameter a is chosen so that 
the frequency 4- vector is K a = dx a / da. Let u denote frequency measured by a stationary 
observer at r. Then we have 

u = U a K a = gooU°K = i^l - 2m/r, ■ = A . 

However (1 — 2m/r)i is constant along the worldline. Therefore wyl — 2m /r is also 
constant and so we have 



6o>2 = U>i > 



1 — 2mjv\ 



1 — 2m /r2 

This is the gravitational redshift. For large r 1; r 2 , we have 

uj 2 ~ o»i(l + m/r 2 - m/ri) , 

so the change in freqeuncy is proportional to the difference gravitational potential between 
the two observers. This is precisely what is needed to avoid the paradox in the working of 
Bondi's perpetual motion machine. 

Remark. The energy of a photon (relative to an observer) is huj. The conservation law 
here can again be interpreted as 'conservation of energy'. 
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12.5 Killing vectors 

A special role is played in these calculations by time symmetry: it is this that allows us 
to say what we mean by 'stationary' observers, and it is this that gives us the energy 
conservation. 

More generally, if the metric coefficients g ab are independent of one of the coordinates 
x°, then L = \g ab x a x b is independent of x°, and so from Lagrange's equations 



is constant along geodesies. But this quantity is equal to T a V a , where V a = dx a /ds and 
T is the 4-vector field with components (1,0,0,0). We note that T a V a is an invariant: it 
depends only on the 4-vectors V and T, and not on the choice of coordinates (although, of 
course, T will have components (1, 0, 0, 0) only for special choices of coordinates). 

Definition 9 

A non-vanishing vector field T is said to be a Killing vector field or Killing vector when- 
ever there exists a coordinate system in which T has components (1, 0, 0, 0) and g ab is 
independent of x°. 

We have just proved the following. 

Proposition 10 

If T is a Killing vector, then T a x a is constant along any geodesic. 

Remark. For any non-vanishing vector field T ^ 0, there exists a coordinate system in 
which T has components (1, 0, 0, 0). 

How can we recognise a Killing vector and therefore derive a conserved quantity for free 
particle and photon orbits without making the transformation to the special coordinate 
system? To answer this, we look first at the defining property in the special coordinates. 
Here we have 



dL 



■a 



dx° 



= gaOX 



= d g a b = T c d c g ab . 



But we also have 



V a T 6 
V„T a 



d a {gb c T c ) - \T c (d a g bc + d b g ac - d c g ab ) 
db{g ac T c ) - \T c {d h g ac + d a g bc - d c g ba ) . 
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By adding, we get 

V a T b + V b T a = T c d c g ab = . 

But the left-hand side is a tensor. Therefore it vanishes in one coordinate system if and 
only if it vanishes in every coordinate system. We have proved the following. 

Proposition 11 

Let T a be a non-vanishing vector field. If T a is a Killing vector then V a T b + V b T a = (in 
any coordinate system). 

The converse is also true: this is proved by choosing a coordinate system in which the 
components are (1, 0, 0, 0), and by using the same calculation in reverse. 

We can use (11) to prove Proposition (10) directly by starting from the geodesic con- 
dition in the form V™V a Vfe = 0, where V a = x a . From this we get 

(V a T a )' = V a V a (V b T b ) = V a V b V a T b = \V a V b (y a T b + V b T a ) = , 

and hence that V a T a is constant. 

Exercise. The converse statement can be deduced from this: if x a T a is conserved along 
every geodesic, then T a is a Killing vector). 

We use Proposition (11) to extend the definition to arbitrary vector fields. We then 
have 

Definition 10 

Alternative defintion: a vector field T a is a Killing vector if V a T b + V b T a = 0. 
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Lecture 13 



13.1 Massive particles 

We shall now look at particle motion in the Schwarzschild background. Our main aim 
will be to derive corrections to Kepler's laws, so we shall think of the gravitational field 
as that of the sun. By a 'particle', we mean a very small body, such a planet, whose own 
gravitational field can be ignored. 

The particle orbits are generated by the Lagrangian 



1-— V 2 ~ 1 r o / -r 2 (v 2 sin 2 "- 2 
r J 1 — 2m I r 



where the parameter is the proper time s and • = d/ds. We shall assume that r > 2m, 
which means that we are looking at the external field of a spherical star rather than the 
field inside a black hole. Since L has no explicit dependence on t, <j), or s, we have three 
conservation laws. 



(d t L = 


0) 


E 


:= (1 — 2m/r)t = constant 


(dfL = 


0) 


J 


:= r 2 sin 2 9 <f> = constant; 


(d s L = 


0) 


L 


= constant. 



Since s is proper time, g a bX a x b = 1, and therefore the third conservation law is L = |. We 
need one other equation to determine the orbits. We shall use the 9 Lagrange equation, 

d ( 'I,-. \ -2 , n ,..„. n '2 



ds V J ~ Sin 6 C ° S 6 ^ ° ' ^ 

We could also write down the r equation, but it would contain no new information since 
we already have four equations for the four unknown coordinates t, r, 9, 0. 

Because (38) is symmetric under 9 i— > tt/2 — 9, an orbit such that 9 = n/2, 9 = 
initially will have 9 = it/ 2 for all s. Because the field is spherically symmetric, we can 
understand all the orbits by studying only these equatorial orbits. So we shall put 9 = tt/2. 
Then by combining the conservation laws, we have 

| E 2 f 2 _ J 2 

1 — 2m I r 1 — 2m jr r 2 
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9 J 2 f 2m\ 9 ( 2m\ 
f 2 = 1 ) + E 1 ' * 



That is, 

r 2 = - s . . . ... 

\ f } V f } 

This is a first-order differential equation for r as a function of proper time. As in Newtonian 
theory, the equation looks a bit simpler if we replace r by u = m/r and use instead of s 
as the parameter. Now 

du mdr / d<f) mr 
d(f) r 2 ds/ ds J 
Hence when J ^ 0, that is when the orbit is not radial, we have 

du\ 2 m 2 E 2 9/ , m 2 (l — 2u) , , 



13.2 Comparison with the Newtonian theory 

In Newton's theory, the particle (assumed to have unit mass) moves under the influence 
of the inverse-square law force m/r 2 . The equatorial orbits are determined in plane polar 
coordinates r, by the conservation of the energy e and the angular momentum J: 

e = \{r 2 + r 2 2 ) — m/r, J = r 2 cf). 

As before, we put u = m/r, du/d(p = —mf/J. Then we have 

_ i / V 

2 I m 2 \d(j) J m 2 J 

To make comparison with Einstein's theory, we put (3 = m/J, p = dw/d0 and k = em 2 / J 2 
(Newton), k = (E 2 — l)m 2 /2J 2 (Einstein). Then we have (from above and from (39): 

Einstein p 2 = 2(3 2 u + 2k-u 2 + 2u 3 =: f(u); 
Newton p 2 = 2/3 2 u + 2k - u 2 := g(u) . 

The only difference is the term 2m 3 , which is, of course, small when r is large. 
The equations for the orbits can also be written in then second-order form: 



d 2 u d ( 1 / du 



2 N 



do* du \ 2 ' ~ 2f ' {u) 

in Einstein's theory, or in the same way with \g'{u) on the right in Newton's theory. 
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13.3 Newtonian orbits and the relativistic correction 



We can see the effect of the extra term in one of the classical tests of general relativity, 
the perihelion advance (i.e. the fact that the point on each orbit where r is minimal — the 
perihelion — advances on each orbit in Einstein's theory, while in the Newtonian theory, it 
is always in the same place. 1 

In the Newtonian theory, we have 

d 2 u m 2 

W + u = ^' (40) 

which implies that u = (3 2 + Acos(0 — 0o), for constant A, </>o. By differentiating we get 

V 2 = A 2 sin 2 (0 - O ) = A 2 - (u - 13 2 ) 2 . 

Hence A 2 = 2k + (3*. 

(1) If k > then \A\ > (3 2 and u = for some values of 0. In this case, the orbits are 
hyperbolic and the particle can escape to infinity. 

(2) If k < then \A\ < (3 2 . In this case, u is bounded away from zero, and therefore |r| is 
bounded and the orbits are elliptic. Special cases are the circular orbits, on which u = uq 
is constant, and so du/d<j) and d 2 u/d<f) 2 vanish identically. These are therefore given by 
solving g(u ) = g'(u ) = for w . The result is u = j3 2 , where 



g(P 2 )=P 4 + 2k = 0. 



That is, k = — 1/3 4 . 

We can rewrite (40) in the form 



d 2 v 



where v = u — (3 2 . This is the equation of simple harmonic motion with period 2tt (the 
'time' of course is not t or s, but the polar angle 0). Thus we can think of a general elliptic 



1 In fact the perihelion of Mercury also advances in Newtonian theory because of interactions with other 
planets, notably Jupiter: the general relativistic effect is the additional advance that cannot be explained 
in this way. When it was first noted, before the discovery of general relativity, it was thought that it might 
be due to another planet orbiting the sun inside Mercury's orbit. This planet, which of course was never 
seen, was called Vulcan. 
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orbit as oscillating about circular orbit (u = (3 2 ) with simple harmonic motion. The fact 
that in these oscillations the period of u as a function of <ft is exactly 2ir is what makes 
the elliptic orbits closed in Newtonian theory: each circuit of the origin adds 2n to and 
brings us back to the initial value of u. In particular, perihelion always occurs at the same 
value of 0: there is no perihelion advance in the two-body system. 

One can gain some insight into the structure of the orbits in Newtonian theory by 
plotting the phase portrait, in which one represents by curves in the w,p-plane. If we fix (3 
and plot the curves for varying values of k, the result is a set of concentric circles 

p 2 + (u - (3 2 ) 2 = 2k + (3 4 

centred on the circular orbit Uq = (3 2 (labelled A). The hyperbolic orbits are those that 
meet the p axis, the elliptic orbits are those that do not. The two families are separated 
by the parabolic orbit, which touches the p-axis at the origin. 
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13.4 The perihelion advance 

In general relativity, there are also closed orbits u = u Q . The corresponding values of the 
constant (3 and E are found from f'(uo) = = f(uo). 

Now consider an orbit u = uq + v(<f>) which is almost circular, so that v is small. By 
substituting into the equation of motion, we obtain 

^ = \ f{u) = \f{u ) + \vj" (u ) + 0(v 2 ) . 

But f'iuo) = 0, and f"(u) = — 2 + 12m. Therefore v satisfies 

+ (1 - 6u )v = 



d0 ; 



(ignoring the 0(v 2 ) error term). This is again the equation of simple harmonic motion 
(with cf) as 'time'); so at least for orbits that are close to circular, we again have the picture 
that the particle's orbit oscillates about a circular orbit. However, not the 'period' is not 
27r, but 

$ = 27T/V1 - Qu ~ 2tt + 6u rc , 

for small «o, that is, large r$. Thus if the particle starts at perihelion (r minimal, u 
maximal), then r returns to its initial value not after a whole rotation, but after (f> has 
advanced through a further angle 6uQn. This is the perihelion advance. 

For small Uq (large r), we have $ ~ 2n + 6^0^; so perihelion occurs at <p = 0,2n + 
6uo7i, . . . (ignoring terms of order u^). If we substitute uq = m/r and put back in the 
constants (there is only one way to do this to get dimensions right), then the advance is 

6Gm7r 



r c 2 



per revolution for an orbit of approximate radius ro (we are ignoring second order terms 
in 1/ro, as we U as assuming that the orbit is 'nearly' circular). 

In the case of Mercury, m, the mass of the sun, is 1.98 x 10 30 . The radius is r = 
5.79 x 10 10 , and the constants are G = 6.67 x KT 11 , and c 2 = 9 x 10 16 (all in SI units). 
This gives the advance as around 40" per century (a more careful analysis gives 43"). In 
the case of the binary pulsar, the advance is around 4° per year. 
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Lecture 14 



14.1 Circular orbits 

The equatorial particle orbits in the Schwarzschild space-time are given by p 2 = f(u), 
where 

f(u) = 2f3 2 u - u 2 + 2u 3 + 2k , 

and u = m/r, p = du/d<p, [3 = m/J, k = \{E 2 — l)m 2 /J 2 . In second order form they are 
solutions to 



2 ' 

d 2 U 



l f{u). 



d<f) 2 2 " 

The circular orbits are those for which r and therefore also u are constant. They are given 
by f( u ) — 0, f'(u) = 0. The second of these implies that 



6u = 1 ± v 7 ! - 12/? 2 , 

which for small (3 gives u = (3 2 and u = | — (3 2 (ignoring terms of order j3 4 ). The first 
root is the Newtonian circular orbit: this is still present in general relativity provided that 
the radius m//3 2 is large in comparison to m. The second is a new feature: it has radius 
close to r = 3m, which is only just above the Schwarzschild radius r = 2m, and of course 
it exists only if the the source of the gravitational field is all contained within the sphere 
r = 3m, so the metric still takes the Schwarzschild form at this radius. We shall see below 
that r = 3m itself is a cicular photon orbit. A particle on the inner circular orbit has to 
be moving close to the velocity of light (relative to a stationary observer). 



14.2 The phase portrait 

We can understand more clearly the pattern of orbits by drawing the phase portrait in the 
w,p-plane for fixed (3 2 and varying k (as we did above for the Newtonian theory). We first 
plot 2f3 2 u — u 2 + 2u 3 for different values of f3 2 : 

(3 2 = 0: (red curve). 

(3 2 = 1/16: the two roots of 2u 2 - u + 2/3 2 come into coincidence (blue curve). 

f3 2 = 1/12: the two roots of the w-derivative of 2f3 2 u — u 2 + 2u 3 come into coincidence 
(green curve). 
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Figure 9: Plots of 2j3 2 u - u 2 + 2u 3 . 



We shall consider the orbits only for | > u > 0, that is, for r > 2m (if the vacuum region 
extnds that far, the portion of space-time in which r < 2m is inside a black hole. 

For different choices of /3 2 , the phase portraits are as follows. Note that for small u 
(large r), the portraits all correspond to the Newtonian one. 

< /? 2 < 1/16. In this case, there are two circular orbits, one stable (5), and the other 
unstable with k > (A). A particle disturbed from the inner cicular orbit — the unstable 
one — can either spiral inwards or escape to infinity. The boundary r = 2m is the red line. 
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1/16 < f3 2 < 1/12. In this case, the inner unstable circular orbit A has k < 0: a particle 
disturbed from this orbit will not escape to infinity. As j3 2 is increased, the two circular 
orbits move towards each other. They coincide when f3 2 reaches 1/12, at r = 6m. 
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f3 2 > 1/12. There are no closed orbits in this case: the angular momentum is too small. 
All orbits either escape to infinity or spiral inwards. 

We note that in no case are there stable circular orbits with r < 6m: this is the minimum 
radius for a planetary orbit. For the sun, this minimum radius is 9 km. That is, it is well 
inside the sun: the limit is not relevant. In the case of black holes, the limit is important in 
the analysis of the infall of matter (usually from a companion star): it is matter falling in, 
and heating up as it does so, that is responsible for X-ray emissions from the neighbourhood 
of stellar mass black hole. 

An interesting lesson to learn from the first case is that, contrary to popular belief, 
it is not easy to fall into a black hole. Suppose that the particle starts a long way from 
the Schrwarzschild radius (r >> m), and that the radial and tranverse components of its 
velocity relative to a stationary oserver are t> rad and v trans , respectively. So long as these 
are small compared with 1 (the velocity of light), we have that 

„ m m 1 
f3=- = 

] ' ^trans 

is small. The corresponding inner circular orbit is given approximately by r = 3m, and so 
2k = 1/27 on this orbit. It is clear from the first phase diagram, therefore, that for the 
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particle orbit to pass inside the black hole (i.e. to pass u = |), it must have 



at large r. But 



p 2 > 



mr 



27 



P 



m v rad 



r 2 r f trans 

If the particle is to fall into the black hole starting at large r, the angle a between its 
trajectory and the inward-point radial direction must satisfy 



^trans SV^fU 

tan a = < 

^rad r 



3v^3 Schwarzschild radius 
2 distance 



So, for example, if the black hole has the mass of the sun (with a Schwarzschild radius of 
3km) and if initially r is the radius of the earth's orbit, then a has to lie in an impossibly 
small range: although the balck hole has a very strong gravitational field, it is very small, 
and is an almost impossible target to hit from any distance. 
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14.3 Photon orbits 



Photon orbits are generated by the same Lagrangian 

but now the dot is differentiation with respect to an affine parameter a. Again we have 
three constants 

Energy E = (1 — 2m/r)i; 

Angular momentum J = r 2 sin 2 9 <f> 

Hamiltonian L = 

only in this case, h = since g a bi a x b vanishes because x a is null. By proceeding in the 
same way as before, we obtain 

p 2 = a 2 + 2u 3 - u 2 , 

where u = m/r, p = dw/d0, and a 2 = m 2 E 2 /J 2 . We now have no sensible Newtonian 
model with which to make comparisons, but we note that without 2m 3 term, the orbits 
would be given by 

u = a cos(0 — 0o), V — ~ a sin(0 — O ) • 

That is, 

r cos(0 — 4>q) = m/a , (41) 

which is the polar equation of a straight line. Thus the term 2u 3 is the gravitational 
contribution: it is responsible for the 'bending of light by gravity'. 

The phase portrait is found in the same way as before. First we plot a 2 + 2m 3 — u 2 
against u for different values of a, as in Figure 14.3 

By taking the square root, we then get the phase portrait in Figure 14.3, with the 
different values of a labelling the different curves in the p, w-space. For large r (small 
u), the orbits look like hyperbolic Newtonian orbits: a photon travelling in from infinity 
will esacpe to infinity, but the trajectory will be deflected. As r decreases, the deflection 
increases, and the orbit can wind around the source many times. Ar r = 3m, it is possible 
for the photon to orbit in an unstable circular orbit. 
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Figure 10: a 2 + 2m 3 — u 2 against u for different values of a 



14.4 The bending of light 

Let us solve for the orbits on the assumption that r is large (and therefore u 2 « u 3 ). 
We shall look at the effect of the 'gravitational term' by putting u = acos<f) + v, where 
a is small and v = 0(a 2 ) (i.e. we are looking at a perturbation of (41) with O = 0). By 
ignoring terms of order v 2 = 0(a 4 ), we get 

= u' 2 + u 2 — 2m 3 - a 2 = -2av' sin + 2av cos - 2a 3 cos 3 + 0(a 4 ) 

where the prime is differentiation with respect to <fi. From this we obtain 



sm 



v = v cos 



— a 2 cos' 



3 




sin sin 
=^> v = a 2 {l + sin 2 0) , 



+ a 2 sin + K 
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Figure 11: The phase portrait for photon orbits 



since without loss of generality, we can set K = by adjusting <j) Q . Since the gravitational 
field is weak, we can interpret r, as plane polar coordinates in the equatorial plane. 
The unperturbed trajectory is the straight line or cos = m: a point on this goes to 
infinity as — > ±7r/2. The perturbed trajectory, on the other hand, goes to infinity as 
— > ±(7r/2 + 7), where the angle 7 is given by 

= —a sin 7 + a 2 (l + cos 2 7) , 

which gives 7 = 2a to the first order in a. The total deflection of the light ray is 27 = 
4a = 4m/ D, where D = m/at is the value of r at the point of closest approach to the 
source of the unperturbed trajectory In SI units, the deflection is 

AmG 
Dc 2 ' 

For a light ray just grazing the surface of the sun, we have (in SI units) D = 7 x 10 8 (the 
radius of the sun), m = 2 x 10 30 (the mass of the sun), c = 3 x 10 8 , G = 7 x 10 -11 . The 
result is a deflection of 10~ 5 radians or 2". This is hard to see, of course, because of the 
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difficulty in observing light that grazes the sun. Its effect can, however, be seen during a 
total eclipse, and was first observed by Eddington in 1918: the deflection causes the stars 
near the sun in sky to move outwards slightly from the centre of the sun, as compared with 
their their normal relative positions. Eddington compared photographs of the star field 
near the sun during a total eclipse with a photograph of the same star field when the sun 
was in a different position in the sky. 
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Lecture 15 



15.1 The Schwarzschild radius 

We shall now look more closely at what happens at the 'Schwarzschild radius', r = 2m. 
Clearly something goes wrong in the formula for the metric coefficients at that point: we 
shall see, however, that it is not the space-time geometry itself that is singular there, but 
simply the coordinates in which it is expressed. We shall see that the singular behaviour at 
r = 2m (but not at r = 0) goes away when we make an appropriate change of coordinates. 

Of course for a normal star the Schwarzschild radius is well inside the star itself: since 
it is not in the vacuum region, we would not expect R a b = to hold at r = 2m, and 
therefore the Schwarzschild solution would not be valid 2 . But, in the extreme case that all 
the body lies within its Schwarzschild radius, we have a (spherical) black hole. 

In the case of the sun, for r = 2m to be outside the sun (and therefore in the vacuum 
region), the sun would have to be compressed to a radius of 3km, which would imply almost 
unimaginable density. For a galaxy, however, the density at this critical compression is only 
that of air, and so it is not hard, at least in principle, to imagine a sufficiently advanced 
civilization directing the orbits of all the stars in galaxy so that all the matter ends up 
within the Schwarzschild radius. We must therefore take seriously the existence of black 
holes theoretical possibility. 



15.2 Eddington— Finkelstein coordinates 

The Schwarzschild metric is 

ds 2 = L _ 2l A dt 2 ^!_^_ r 2 (d0 2 + sin 2 0#2) 

\ r J 1 — 2m/r 

We cannot simply ignore the part of space-time for which r < 2m because an infalling 
observer will reach r = 2m in finite proper time. The observer's worldline (if he falls 
radially, that is with constant 9 and 0) is given by 

E = (1 - 2m/r)t, 1 = (1 - 2m/r)i 2 - -, 

1 — 2m jr 

2 One can construct 'interior' Scwharzschild metrics with the energy-momentum tensor of various types 
of matter on the right-hand side of Einstein's equation; these are, of course, nonsingular at r = 2m. 
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where the parameter s is proper time. If we consider the special case E — 1, which arises 
when the observer falls from rest (with respect to the time Killing vector) at infinity, then 
we have r 2 = 2m/ r. For an ingoing observer, therefore, 



rdr = — v2m J ds 

_^ | r3//2 = \ / 2m(k — s) ; k = constant . 

from which we conclude that the proper time s taken to reach r = 2m is finite. However 
the coordinate time taken, t, is infinite because 

dr / 2m\ \/2m 



dt \ r J \Jr 

and so 

/2m I dt. 



r — 2m 

But the intergal on the left-hand side diverges as r — > 2m. 

To understand the space-time geometry of a black hole, we first look for a coordinate 
system in which the singularity at r = 2m disappears. One can see what goes wrong with 
the given coordinates by looking at the null geodesies in the r, t-plane (the worldlines of 
photons travelling radially inwards or outwards). These are the curves given by 

2m\ , 9 dr 2 



1 dt 2 = 

r J 1 — 2m /r 



from which we obtain 



r dr f ( 2m \ 

dt = ± = ± / 1 + dr . 

J 1 — 2m / r j \ r — 2m J 

That is, 

t ± (r + 2m log(r — 2m)) = constant . (42) 

The radial null geodesies in the t,r-plane are shown in blue in Figure 15.2. The singular 
behaviour at r = 2m (the red line) is associated with the fact that the curves 'bunch up' 
on the common asymptote at r = 2m. We also see that for large r, they look like the 
corresponding curves r = ±t in flat space-time. One should think of each curve as an 
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Figure 12: Radial null geodesies in the Schwarzschild metric 

ingoing or outgoing spherical wave front: one can get a partial picture of this by rotating 
the about the t-axis to make the curves into surfaces of revolution. The are shown Figure 
15.2, which is a space-time diagram with one spatial dimension suppressed. The blue 
surface is at r = 2m; the history of the outgoing wave-front is shown in red; and that of 
the ingoing one in green. 
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Figure 13: Ingoing and outgoing wavefronts in the Schwarzschild metric 
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Figure 14: Radial null geodesies in Minkowski space 

We should compare these pictures with the corresponding ones for Minkowski space. 
The corresponding radial null geodesies are the straight blue lines shown in Figure 15.2, 
and the corresponding in- and outgoing wave-fronts are the null cones of the points on the 
polar axis r = (Figure 15.2). 
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Figure 15: Ingoing and outgoing wavefronts in Minkowski space 



In the Schwarzschild space-time, we can resolve the coordinate difficulties by 'compress- 
ing' the t coordinate as we approach r = 2m. Guided by (42), we make the transformation 
to coordinates v, r, 9, by putting 

v = t + r + 2m log(r — 2m) , 

which gives 

dr 

dt = dv 



1 — 2m jr 
and hence 

ds 2 = (1 - 2m/r) dv 2 - 2dv dr - r 2 (d6 2 + sin 2 9 d0 2 ) . 

The singular behaviour at r = 2m has now disappeared. In the r, f-plane, the radial null 
geodesies are given by 

2m \ dv 

v = constant or 1 2 = 0. 



r J dr 
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The second equation can be integrated to give: 



v — I = 2r + 4m log |r — 2m\ + const. 

I r — 2m 



Thus the radial null geodesies are as shown in Figure 15.2. 




Figure 16: Radial null geodesies in Eddington-Finkelstein coordinates 

We see from the from the fact that timelike curves must lie between the ingoing and 
outgoing null geodesies at every event that, although the space-time is nonsingular for 
r < 2m, it is not possible to escape to infinity. Hence the hypersurface r = 2m is called the 
'event horizon': it separates events of which observers outside can have knowledge from 
those inside of which they cannot. The events inside the event horizon are inside the 'black 
hole'. (The v-axis has been drawn at 45° to emphasize the fact that the lines of constant 
v are null.) 

The histories of the ingoing and outgoing wave fronts outside the event horizon(the 
blue surface) are now as in Figure 15.2. 

The singular behaviour of the the metric coefficients in the t, r coordinates is not a 
singularity of the space-time geometry (because it disappears in the v,r coordinates): it 
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Figure 17: Ingoing and outgoing wavefronts in Eddington-Finkelstein coordinates 

arises from the singular behaviour of the transformation from v,r to t, r coordinates at 
r = 2m — which shows itself in the picture in the fact that the curves of constant t are 
asymptotic to the line r = 2m. The transformation from v, r to t, r coordinates pushes the 
points of (v, 2m) to t = oo. 

The Killing vector T with components (1,0,0,0) in the coordinates t,r,6,<p has the 
same components in the new coordinate system, but we note that T a T a = 1 — 2m jr both 
inside and outside the Schwarzschild radius. However 

(i) for r > 2m, T is timelike and defines a standard of 'rest' (an observer whose 
4- velocity is tangent to T is what we called a 'stationary observer'); 

(ii) for r = 2m, T is null: we can think of the event horizon as the history of 
a light wavefront 'at rest', hovering for ever between escaping to infinity and 
falling into the black hole; 

(iii) for r < 2m, T is spacelike, and no observer can remain at rest. 

The worldline of any observer inside the black hole must inevitably reach r = in finite 
proper time (in fact, in a time of the same order of magnitude as light takes to travel a 
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distance equal to the Schwarzschild radius). In fact r = is a genuine singularity at which 
the tidal forces become infinite. One can see this from the fact the invariant R a b c dR abcd 
blows up like r~ 6 , and therefore there cannot exist a coordinate system in which the metric 
is well-behaved at r = 0. Thus once inside the black hole, an observer is not only unable 
to escape to infinity but is equally unable to escape being crushed in the singularity in 
a very short time (the position is rather different if the black is rotating: see next term's 
lectures) . 
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Lecture 16 



16.1 Gravitational collapse 

The Schwarzschild solution by itself does not provide a good model of a real black hole 
because it is a vacuum metric: there is no matter present to generate the gravitational 
field. In a real astrophysical situation one expects black holes to form from the collapse 
of stars after they have burnt up all their nuclear fuel: the collapse can form a white 
dwarf, which is supported against gravity by the 'electron degeneracy pressure'; however, 
above 1.4 times the mass of the sun, this pressure is insufficient, and collapse results in 
a neutron star — essentially a massive nucleus with an atomic number around 10 58 . But 
again there is a limit to mass: above some critical mass (somewhere between 1.5 and 3 
solar masses), no known physical process can prevent collapse to a black hole. Once the 
event horizon has formed, no conceivable process can prevent collapse to a singularity (the 
Penrose singularity theorem). 

One can model the field of a spherically symmetric collapsing object by joining the 
Schwarzshild metric (to represent the field outside the body) to an interior metric (rep- 
resenting the field inside the collapsing star) across a spherically symmetric hypersurface 
represented by a timelike curve in the w,r-plane. 

If we add one of the other coordinates by rotating about the line r = 0, then one obtains 
the following three-dimensional representation of the space-time: 

16.2 Kruskal coordinates 

It is interesting, however, to explore further the vacuum solution without joining on any in- 
terior solution. Here we look more closely at a curious feature of the Eddington-Finkelstein 
coordinates: they introduce time asymmetry that is not present in the original metric. 
That is, they do not treat the future and the past in an even-handed way. They adjoin 
the interior of a black hole to the exterior solution, or in the time-reversed form, a 'white 
hole'. 

One can see what is going on here by transforming instead to Kruskal coordinates, in 
which both extensions can be made simultaneously. We start with the original form of the 
metric 

ds 2 = ( 1 _ 2 lA dt 2 ^^_ r 2 (d 02 + gin 2 #2) 

\ r J 1 — 2m/r 
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Figure 18: The collapse of a star to form a black hole 

But now we transform to new coordinates U, V, 6, <p by putting 

V 



u 



-e t/2m , UV = e r/2m (2m - r) 



That is, V = e"/ 4m , U = - e ~ u ^ m , where 

v = t + r + 2m log(r — 2m), u = t — r — 2m log(r — 2m) 
(Note that v is the Eddington-Finkelstein coordinate). Then 

gf/4m / dr \ Q—u/im / ^ r 



dV = ~ I dt , 

4m \ 1 — 2m I r 



dU 



Hence 



16m z 



r 



4m \ 1 — 2m /r 
dr 2 



1 — 2m/r) 



Therefore in these new coordinates, the metric is 



ds 2 = 16m 2 r- 1 e- r/2m dEW - r 2 (dfl 2 + sin 2 9 d(j> 2 ) 
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where r is defined as a function of U, V by UV = e r/2m (2m - r). 

To understand the geometry, let us look first at the corresponding transformation of 
Minkowski space. Here we start with 

ds 2 = dt 2 - dr 2 - r 2 (d# 2 + sin 2 6 d0 2 ) 

and make the coordinate change 

U = -e r ~\ V = e t+r . 

Then the metric becomes 

ds 2 = e- 2r d£/dV-r 2 (dfl 2 + sin 2 fld0 2 ) 
= -^^-r 2 (d£ 2 + sin 2 £d0 2 ). 

If we suppress the angular coordinates, then the relationship between the two coordinate 
systems as shown in Figure 16.2. The U, V axes are null lines, and are therefore drawn at 




Figure 19: The U, V coordinates on Minkowski space 
45° to the horizontal, with time and the two coordinates U, V increasing up the page. The 
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curves of constant t (green) are straight lines through the orgin; those of constant r are the 
hyperbolas UV = constant, which have the U, V -clXCS clS asymptotes. The transformation 
maps the whole of Minkowski space into the region 

-UV > 1, U<0, V>0 

in the U, V-pla.ne. The blue hyperbola is the curve UV = —1; that is, r = 0. 

In the Schwarzschild geometry, the picture is very similar, except that the metric con- 
tinues in the U, V-pl&ne to the region UV < 2m (the boundary UV = 2m is the image of 
the 'real' singularity at r = in the r,t plane). Figure 16.2 again shows the U, V^-plane, 




Figure 20: The Krusakal extension of the Schwarzschild geometry 

with the axes drawn at 45° to the horizontal. The blue lines are null. In this case, how- 
ever, the metric is non-singular in the whole region bounded by the two branches of the 
hyperbola UV = 1 (on which r = 0): if we exclude the red-shaded region above and below 
these, then we have the maximal analytic extension of the Schwarzschild space-time, due 
to Kruskal. 

If the black hole is formed by gravitational collapse, then we see only part of the diagram 
to the right: the rest must be replaced by a suitable interior metric. The portion covered by 
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the Eddington-Finkelstein coordinates is the portion above the [/-axis. The entire extended 
space-time contains both a black hole (the region U > 0, V > 0, which an observer can 
enter but not leave), and a 'white hole' — the time reverse of a black hole — namely the 
region U < V < 0, which an observer can leave but not enter. 

If we admit all values of U, V with UV < 2m, then we have the 'Kruskal extension' of 
the Schwarzschild space-time, which contains no matter: we can think of the 'm' in the 
metric as being entirely gravitational in origin, or perhaps we should think of it as the 
mass of the singularity at r = 0. Here we have no stellar boundary, and space-time looks 
like two 'external' regions, but they are joined together by a 'wormhole'. 

The external regions are the two quadrants V > > U and U > > V: for large \UV\, 
the metric looks in both like that of Minkowski space . We can see the way in they are 
connected by looking at the geometry of the spatial section U — V, on which r is given as 
a function of V by V 2 = e r / 2m (2m — r). 




Figure 21: The spatial geometry at t — 

On this r decreases to a minimum value of 2m and then increases again to infinity. If 
we put 9 = 7r/2 (so that we are looking at the 'equatorial plane'), then the metric is 

ds 2 = (1 - 2m/r)- 1 dr 2 + r 2 d0 2 = (1 + f(r) 2 )dr 2 + r 2 d0 2 , 



86 



where / = ^/8m(r — 2m). This is the metric on a surface of revolution given by rotating 
the parabola / = f(r) about the /-axis (Figure 16.2). Thus we can picture the hypersurface 
U = V as two copies of Euclidean space (at large r), joined by a wormhole: to an observer 
in either of these external spaces, the geometry looks like that of a black hole. Of course 
one cannot actually travel through the wormhole: passing through r = 2m takes one inside 
the event horizon, and inevitably into the singularity at r = 0. 



87 



